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ABSTRACT: Observations of damped Lya systems offer a unique window on the neutral-gas 
reservoirs that gave rise to galaxies at high redshifts. This review focuses on critical properties 
such as the H I and metal content of the gas and on independent evidence for star formation. 
Together, these provide an emerging picture of gravitationally bound objects in which accretion 
of gas from the IGM replenishes gas consumed by star formation. Other properties such as dust 
content, molecular content, ionized-gas content, gas kinematics, and galaxy identifications are 
also reviewed. These properties point to a multi-phase ISM in which radiative and hydrodynamic 
feedback processes are present. Numerical simulations and other types of models used to describe 
damped Lya systems within the context of galaxy formation are also discussed. 
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1 WHAT ARE THE DAMPED Lya SYSTEMS and HOW ARE 
THEY FOUND? 

Damped Lya systems are a class of QSO absorbers selected for the presence of H I 
column densities, N(H I)> 2xl0 20 cm~ 2 . This criterion differs from those used to 
find other classes of QSO absorbers selected on the basis of H I content. The Lya 
forest absorbers, reviewed in this journal by Rauch (1998), are selected for N(H I) 
< 10 17 cm~ 2 , while the Lyman limit systems have 10 17 < iV(H I) < 2xl0 20 cm -2 
(Peroux et al. 2003b). The Lya forest absorbers are optically thin at the Lyman 
limit, since the column density N(R I) = 10 17 cm' 2 corresponds to about unity 
optical depth at the Lyman limit. Are these absorbers physically different from 
the damped systems or have the column-density criteria resulted in arbitrary 
distinctions? In fact there is a fundamental difference: hydrogen is mainly neutral 
in damped Lya systems, while it is ionized in all other classes of QSO absorption 
systems. This includes absorbers selected for the presence of C IV AA 1548.1, 
1550.7 resonance-line doublets (Sargent, Steidel, & Boksenberg 1988), Mg II AA 
2796.3, 2803.5 resonance-line doublets (Steidel & Sargent 1992) and Lyman limit 
absorption (Prochaska 1999), which do not also qualify as damped Lya systems. 

The neutrality of the gas is crucial: while stars are unlikely to form out of warm 
ionized gas, they are likely to descend from cold neutral clouds, which are the 
precursors of molecular clouds, the birthplace of stars (Wolfire et al. 2003). This 
property takes on added significance when it is realized that the damped Lya 
systems dominate the neutral-gas content of the Universe in the redshift interval 
z=[0,5], and at z~3. 0—4.5 contain sufficient mass in neutral gas to account for 
a significant fraction of the visible stellar mass in modern galaxies (e.g. Storrie- 
Lombardi & Wolfe 2000). This has led to the widely accepted idea that damped 
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Lya systems serve as important neutral gas reservoirs for star formation at high 
redshifts (e.g. Nagamine, Springel, & Hernquist 2004a). Moreover, as repositories 
of significant amounts of metals the damped Lya systems have been used to trace 
the age-metallicity relationship and other aspects of galactic chemical evolution 
(Pettini et al. 1994; Pei, Fall, & Hauser 1999; Pettini 2004; Prochaska et al. 
2003a). 

The purpose of this review is to present an overview of the damped Lya sys- 
tems. Current research on the high-redshift Universe is dominated by surveys 
that rely on the detection of radiation emitted by stars (e.g. Steidel et al. 2003; 
Giavalisco et al. 2004; Dickinson et al. 2003) or ionized gas (e.g. Rhoads & Malho- 
tra 2001; Ouchi et al. 2003). By contrast, damped Lya systems provide a window 
on the interplay between neutral gas and newly formed stars, i.e., the damped 
Lya systems are the best, perhaps the only, examples we have of an interstellar 
medium in the high-redshift Universe. Consequently, the focus of this review will 
be on the manner in which damped Lya systems trace, and play an active role 
in, cosmic star formation and hence galaxy formation. 

Throughout this article we adopt a cosmology consistent with the WMAP 
(Bennett et al. 2003) results, (O m , tt A ,h) =(0.3,0.7,0.7). 

1.1 History of damped Lya Surveys 

To understand the significance of damped Lya systems for research in galaxy 
formation we give a brief historical perspective. 

The motivation for the first damped Lya survey was to find the neutral-gas 
disks of galaxies at high redshifts (Wolfe et al. 1986). Unlike today, the cold 
dark matter (hereafter CDM) paradigm of hierarchical structure formation (i.e., 
merging protogalactic clumps) did not dominate theories of galaxy formation in 
the early 1980s. Rather, the idea of mature galaxy disks at high redshift fitted 
in with the coherent collapse model of Eggen, Lynden-Bell, & Sandage (1962, 
also Fall & Efstathiou 1980), which was highly influential at the time. Some 
QSO absorbers with properties resembling galaxy disks had been found at z < 
1 through the detection of 21 cm absorption either in radio-frequency scanning 
surveys (Brown & Roberts 1973) or at the redshifts of Mg II selected absorbers 
(Roberts et al. 1976). However, application of these techniques resulted in only 
a few detections. While it was unclear whether the 21 cm absorbers belonged to 
a new population of objects or were rarely occurring oddities, the radio scanning 
techniques were valuable for successfully detecting cold, quiescent gas at large 
redshifts for the first time. Specifically, Brown & Roberts (1973) and Brown 
&; Mitchell (1983) used this technique to detect two 21 cm lines with FWHM 
velocity widths, Avui ~ 10 and 20 km s _1 . The temperature of gas detected 
in 21 cm absorption is likely to be low, because the 21 cm optical depth T21 oc 
N(H I)/(T s Avhi), where the hyperfine spin temperature, T s , generally equals the 
kinetic temperature of the cold, dense gas detected in 21 cm absorption. 

However, the most efficient method for locating quiescent layers of neutral gas 
is through the detection of damped Lya absorption lines. In the rest frame of 
the atom the absorption profile of any atomic transition is naturally broadened 
owing to the finite lifetime of the upper energy state. In the rest frame defined by 
the average velocity of the gas, the natural profile is Doppler broadened by the 
random motions of the atoms: the convolution of both effects results in the Voigt 
profile (e.g. Mihalas 1978). Because the Doppler profile falls off from the central 
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frequency, vq, as exp[— {Av/ Avjj) 2 ] (where Au = \v— vq\ and Aur>=y/2a v u /c for 
an assumed Gaussian velocity distribution with dispersion o v ) and the natural 
or "damped" absorption profile falls off from v Q like l/(Av) 2 , at sufficiently large 
Au the probability for damped absorption exceeds the probability for absorp- 
tion in the Doppler profile. The frequency intervals in which natural broaden- 
ing dominates Doppler broadening are called the damping wings of the profile 
function. Most atomic transitions of abundant ions are optically thin in their 
damping wings but optically thick near the core of the Doppler profile. The lat- 
ter transitions have unit optical depth at Av T= \ oc Ai/px [In NiX?)] 1 / 2 , where 
iV(X : ') is the column density of ionic species X- 7 . Such lines are saturated. The 
reason is that the rest-frame equivalent width of an absorption line is given by 
W r =(X/v) J(l — exp(—T u )dv), and therefore W r is proportional to Au T= ±. In the 
case of lines with unit optical depth near the Doppler core the line is saturated 
because W r is insensitive to the value of iV^X-?). Due to the higher values of N(H 
I), Lya has unit optical depth in the damping wings at Av T= \ oc L4.2i/2i-/V(H 
I)} 1 / 2 when N{R I) > 10 19 cm" 2 and u v < 70 km s" 1 : A 2 i and f 2 i are the Einstein 
spontaneous emission coefficient and oscillator strength for the Lya transition. 
In this case, unit optical depth occurs in the damping wings, and therefore the 
equivalent width of a damped Lya line is independent of the velocity structure of 
the gas for velocity dispersions within the range detected in most QSO absorption 
systems. As a result, the equivalent width will be large even when the velocity 
dispersion is small. 

By the early 1980s only four damped Lya systems had been found. In every 
case they were high column-density systems, N(R I) > 10 21 cm~ 2 , which were 
found by chance (Beaver et al. 1972; Carswell et al. 1975; Smith, Margon, & 
Jura 1979; Wright et al. 1979). Although the sample was sparse, the utility 
of the damped Lya criterion was demonstrated when 21 cm absorption at z~2 
was detected in 2 of the 3 background QSOs that were radio sources (Wolfe & 
Davis 1979; Wolfe et al. 1981). The narrow line widths, Ai>hi~20 km s _1 (where 
AfHi=v / 81n2<7„), and relatively low spin temperatures, T s < 1000 K, implied that 
these absorbers were H I layers in which the gas was cold and quiescent. 

For these reasons Wolfe et al. (1986) began a survey for damped Lya systems 
by acquiring spectra of large numbers of QSOs and then searching them for the 
presence of damped Lya absorption lines. The survey for damped Lya systems 
had several advantages over surveys for 21 cm absorption lines. For example the 
redshift interval covered by a single optical spectrum, Az ~ 1, is large compared to 
that sampled by bandpasses then available for 21 cm surveys, Az 0.02. Second, 
optical spectra of QSOs are obtained toward continuum sources with diameters 
less than 1 pc, whereas the diameters of the associated background radio sources 
typically exceed 100 pc at the low frequencies of redshifted 21 cm lines. As a result 
the survey was capable of detecting compact gaseous configurations with low 
surface covering factors that would have been missed in 21 cm surveys. Another 
advantage of optical surveys is the large oscillator strength, /21 =0.418, of the Lya 
transition (by comparison /21 = 2.5xl0 -8 times a stimulated emission correction 
of 0.068K/T S for the 21 cm line), which allows for the detection of warm H 
I, which is optically thin to 21 cm absorption owing to high values of T s but 
optically thick in Lya. But this is also a disadvantage: the strength of the Lya 
transition combined with the high abundance of hydrogen means that the more 
frequently occurring low column-density clouds in which H is mainly ionized will 
be optically thick in Lya. The result is a profusion of Lya absorption lines, 
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i.e., the Lya forest, which dominate the absorption spectrum blueward of Lya 
emission (see Figure 1). Although the Lya forest lines act as excellent probes 
of the power spectrum and other cosmological quantities (see Tytler et al. 2004; 
McDonald 2003), they are potential sources of confusion noise for the detection 
of damped Lya lines, especially at z > 4, since the line density per unit redshift 
increases with redshift. Identification of damped Lya lines at z > 5.5 is essentially 
impossible because of Lya forest confusion noise. 

However, at z < 5.5 the large column densities of H I in galaxy disks or in any 
other configuration produce damped Lya absorption lines that are strong enough 
to be distinguished from the Lya forest (see Figure 1). Consider the equivalent 
widths. At the time of the Wolfe et al. (1986) survey the most accurate 21 cm 
maps of spiral galaxies were obtained with the Westerbork radio interferometer. 
These showed the H I column densities of galaxy disks to decrease from iV(H 
I) ~ 10 21 cm~ 2 at their centers to iV(H I)= 2xl0 20 cm -2 at a limiting radius 
i?;=(1.5±0.5)i?26.5) which was set by the sensitivity available with Westerbork 
and comparable radio antennas. Here the Holmberg radius, R26.5, is the radius 
at which the B band surface brightness equals 26.5 mag arcsec -2 (Bosma 1981). 
The rest-frame equivalent width of a damped Lya line created by an H I column 
density, N(H I), is given by W r « 10x[7V(H I)/2xl0 20 cm~ 2 ] 1 / 2 A. Because the ob- 
served equivalent width of a line formed at redshift z is W \ )S ={l+z)W r , damped 
Lya systems with N(H I)>2xl0 20 cm will appear in optical QSO spectra with 
Wobs > 16A for damped Lya systems redshifted redward of the atmospheric cut- 
off (i.e. z > 1.6 for X atm = 3200A). Lines this strong are easily distinguishable 
from the W b. s ~ 3 A equivalent widths of typical Lya forest lines. Further- 
more, they can be detected at low resolution and moderate signal-to-noise ratio. 
Since the goal of the first survey for damped Lya systems was to find absorbers 
with A^(H I) > 2xl0 20 cm -2 , a spectral resolution, AA =10 A, was sufficient for 
resolving candidate features. 

1.2 Modern Surveys and Identification of damped Lya systems 

Since the initial survey was published, nine more surveys have been completed for 
damped Lya systems with iV(H I) > 2xl0 20 cm~ 2 (Lanzetta et al. 1991; Lanzetta, 
Wolfe, & Turnshek 1995; Wolfe et al. 1995; Storrie-Lombardi & Wolfe 2000; Rao 
& Turnshek 2000; Ellison et al. 2001; Peroux et al. 2003b; Prochaska & Herbert- 
Fort 2004;Prochaska, Herbert-Fort, & Wolfe 2005). The identification of damped 
Lya systems is more complex than for other classes of QSO absorbers. The 
Lya forest " clouds" , which dominate the absorption spectrum blueward of Lya 
emission, are abundant and easy to identify. Similarly, surveys for C IV or Mg 
II absorption systems rely on the detection of doublets with known wavelength 
ratios, which are straightforward to locate redward of Lya emission. By contrast, 
the task of surveys for damped Lya systems is to pick out a single damped Lya 
line from the confusion noise generated by the Lya forest. In particular, one 
must distinguish a single, strong Lya absorption line created in high column- 
density gas with low velocity dispersion, but broadened by radiation damping, 
from strong Lya absorption features that are Doppler-broadened blends of several 
lines arising from redshift systems with low column-density gas. The presence of 
narrow Lya forest absorption lines in the damping wings of the absorption profile 
is a further complication which can distort the shape of the true line profile in 
data of moderate or low signal-to- noise ratios (see Figure 2 for examples) . 
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The most widely used strategy for discovering damped Lya systems was first 
introduced by Wolfe et al. (1986) and later refined by Lanzetta et al. (1991) and 
Wolfe et al. (1995). First, a continuum is fitted to the entire QSO spectrum 
blueward of Lya emission. Then damped Lya candidates are identified as ab- 
sorption features with rest equivalent widths W r exceeding Wthresh = 5A . This 
conservative criterion corresponds to iV(HI) > 5 x 10 19 cm~ 2 , which guarantees 
that few systems with iV(H I) above the completeness limit of 2xl0 20 cm~ 2 will 
be missed. The search is carried out in the redshift interval z=[z m i n , z max ] where 
Zmin generally corresponds to the shortest wavelength for which a(W r ) < I A and 
Zmax is set 3000 km s _1 below z em to avoid contamination by the background 
QSO. Finally, a Voigt profile is fitted to the Lya profile to determine the value of 
N(K I). Where possible, the centroid is identified from the redshift determined 
by metal lines outside the Lya forest. This is particularly important at z > 3 
where line-blending from the Lya forest often contaminates the damping wings 
(e.g. Figure 2). The surveys were time consuming because the signal-to-noise 
and resolution of the spectra used to acquire damped Lya system candidates 
were usually inadequate for fitting Voigt profiles to the data. Therefore, follow- 
up spectroscopy at higher spectral resolution and with longer integration times 
was usually necessary. 

Recently Prochaska & Herbert-Fort (2004) and Prochaska, Herbert-Fort, & 
Wolfe (2005) have streamlined this process in a survey based on a single set 
of QSO spectra drawn from the Sloan SDSS archive (Abazajian et al. 2003). 
Because of the high quality, good spectral resolution (R ~ 2000) and extended 
spectral coverage of the data, the authors could fit accurate Voigt profiles to the 
same data used to find damped Lya system candidates. The authors also bypass 
the time-consuming step of fitting a continuum to the QSO spectrum blueward 
of Lya emission by searching for damped Lya system candidates in spectral 
regions with lower-than-average signal-to-noise ratios; i.e., regions coinciding with 
broad absorption troughs. The survey is not formally complete to JV(H I) = 
2xl0 20 cm~ 2 , but the similarity between dAf/dX, the number of damped Lya 
systems encountered per unit absorption distance along the line of sight (see 
§ 2.1), in their survey and previous surveys suggests they are more than 95% 
complete. The number of damped Lya systems for the SDSS DR2 and DR3 
archives is 525. As a result the number of damped Lya systems in a statistically 
complete sample now exceeds previous samples by an order of magnitude at z ~3 
and several times at z ~4. 

While the H I selection methods are successful at finding damped Lya systems 
at z > 1.6, they have been unsuccessful at finding large numbers of objects 
at lower redshifts. This is partly due to the reduced interception probability 
per unit redshift at low z and partly because few QSOs have been observed 
from space at UV wavelengths, which is required to detect Lya at z < 1.6. To 
increase the number of low-redshift damped Lya systems from the two confirmed 
objects detected in previous H I selected surveys (see Lanzetta et al. 1995), Rao & 
Turnshek (2000) searched for damped Lya systems in samples of QSO absorption 
systems selected for Mg II AA 2796.3, 2803.5 absorption. Since Mg II absorption 
is present in every damped system in which it could be observed, it turns out 
to be a reliable indicator for the presence of damped Lya. Using this technique, 
Rao, Turnshek, and collaborators have recently increased the sample size to 41 
damped Lya systems with z < 1.6 (SM Rao, DA Turnshek, & DB Nestor 2004, 
priv. comm.). 
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The current sample of damped Lya systems that are drawn from surveys with 
statistically complete selection criteria comprises over 600 redshift systems. While 
the number of damped Lya systems is smaller than the ~ 2350 objects comprising 
the population of known Lyman Break galaxies (Steidel et al. 2003), we expect 
the damped Lya system population to approach this number when all QSO 
spectra from the Sloan database become available. 

1.3 The Significance of the iV(H I) > 2xl0 20 cm~ 2 Survey Threshold 

The survey statistics cited above refer only to systems with N(H I) > 2xl0 20 cm~ 2 , 
which is an historical threshold set by the H I properties of nearby spiral galax- 
ies (see § 1.1). Because the nature of damped Lya systems is still not under- 
stood, their H I properties may differ from those of nearby H I disks: for exam- 
ple, CDM cosmogonies envisage damped Lya systems as merging protogalactic 
clumps (Haehnelt, Steinmetz, & Rauch 1998). As a result, it is reasonable to ask 
whether the 2xl0 20 cm -2 threshold is the appropriate one. Indeed, since the em- 
pirically determined frequency distribution of H I column densities increases with 
decreasing iV(H I) (see Figure 3), lower H I thresholds would be advantageous 
because they would result in larger samples. 

Fortuitously, the 2xl0 20 cm threshold is optimal for physical reasons unre- 
lated to the properties of galaxy disks. Rather, at large redshifts this is the column 
density which divides neutral gas from ionized gas: at iV(H I) < 2xl0 20 cm -2 the 
gas is likely to be ionized while at iV(H I) > 2x 10 20 cm~ 2 it is likely to be neutral. 
The minimal source of ionization is background radiation due to the integrated 
population of QSOs and galaxies. Using background intensities computed by 
Haardt k Madau (1996, 2003), Viegas (1995) and Prochaska & Wolfe (1996) 
show that the gas in most of the "sub-damped Lya" population (defined to have 
10 19 <iV(H I) < 2xl0 20 cm- 2 ) described by Peroux et al. (2002, 2003a) is in fact 
significantly ionized with temperature, T > 10 4 K. This is a problem since gas 
neutrality is a necessary condition if damped Lya systems are to serve as neutral 
gas reservoirs for star formation at high redshift, a defining property of the pop- 
ulation. For this reason the comoving density of H I comprising the sub-damped 
Lya population discussed by Peroux et al. (2003b) should not be included in 
the census of gas available for star formation. As a result, the sub-damped Lya 
correction to the comoving density of neutral gas, Q g (z), should be ignored. We 
suggest that these ionization levels make "super Lyman-limit system" a more 
appropriate name for systems with 10 19 < iV(H I) < 2xl0 20 cm~ 2 . 

2 THE NEUTRAL GAS CONTENT OF THE UNIVERSE 

In this section we describe how the surveys allow us to measure Q g (z), the mass 
per unit comoving volume of neutral gas in damped Lya systems at redshift z 
divided by the critical density, p C rit- The results, first derived by Wolfe (1986) 
and Lanzetta, Wolfe, & Turnshek (1995) show that damped Lya systems contain 
most of the neutral gas in the Universe at redshifts 1.6 < z < 5.0. 

2. 1 Formalism 

To estimate £l g {z) we first derive an expression for the column density distri- 
bution, f(N,X). Let the number of absorbers per sightline with H I column 
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densities and redshifts in the intervals (N, N + dN) and (z, z + dz) be given by 

dAf{N,z) =n co {N,z)A(N,z){l + zf\cdt/dz\dNdz , (1) 

where n co (N, z)dN is the comoving density of absorbers within (iV, TV + dN) at 
z and A(N, z) is the absorption cross section at (N, z). Defining dX=(Ho/c)(l + 
z) 3 \cdt/dz\dz (Bahcall & Peebles 1969) we have 

m ^-! N ~mi(x,N ) , (2) 



dX 



where 



f(N,X)=(c/H )n co (N,X)A(N,X) , (3) 

and where N m i n and N max are minimum and maximum column densities 1 . There- 
fore, one cannot infer the comoving density nor the area of damped Lya systems 
from their incidence along the line of sight, but only their product. Note, dN jdX 
will be independent of redshift if the product of the comoving density and absorp- 
tion cross section at (N, X) is independent of redshift. Since the gaseous mass 
per damped Lya system is given by pmnN A(N : X), it follows from eq. (3) that 

n floMmn f N - dNNmx) . (4) 

C Pcrit JN min 

where fx is the mean molecular weight, which is included to account for the 
contribution of He to the neutral gas content. 

Using these expressions in the discrete limit, several authors have determined 
f(N, X) and its first two moments, dN jdX and £l g (z), where 



Q g ( X ) = ^_ vu , (5) 



Hp pm H ELi N j 

C p cr it AX 

and n is the number of damped Lya systems within (X, X + AX). We now 
discuss each of these in turn. 



2.2 f(N,X) 

Figure 3 shows the most recent determination of f(N,X) from the statistical 
sample of over 600 damped Lya systems (Prochaska, Herbert-Fort, & Wolfe 
2005). The figure also shows best-fit solutions for the three functional forms 
used to describe f(N,X): a single power-law, f(N,X) =k\N ai ; a T function 
(e.g. Pei & Fall 1993) f(N, X)=k 2 (N/N 1 ) a2 exp(-N/N 1 ); and a double power- 
law f(N,X)=k 3 (N/N 1 ) (3 where (3= a 3 at N < N d and a 4 at N > N d . The 
single power-law solution with a best-fit slope of oi\ = — 2.20±0.05 is a poor 
description of the data since a KS test shows there is a 0.1% probability that the 
data and power-law solution are drawn from the same parent population. This 
result is in contrast with the Lya forest where a single power-law with a.\ ~ —1.5 
provides a good fit to the data (Kirkman & Tytler 1997). 

Although a single power-law is a poor fit to the observations, the f(N,X) 
distribution is steeper than iV~ 2 at large column densities. This is illustrated 
by the other two curves in Figure 3 that show the T function (dashed line) and 

^ote that dX/dz = (1 + zf[(l + z) 2 (l + zfl m ) - z(z + 2)n A ] _1/2 
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the double power-law (dashed-dot line). Both solutions are good fits to the data. 
Furthermore, the solutions provide good agreement between the 'break' column 
densities 7V 7 and Nj, and between the power-law indices at low column densities, 
which approach a 'low-end' slope a=— 2.0. Most importantly, both solutions 
indicate a << —2.0 at N > 10 21 ' 5 cm" 2 : the significance of this very steep 
slope at the 'high end' will be explored further in § 2.4. 

Prochaska, Herbert-Fort, & Wolfe (2005) also find evidence for evolution in 
f(N,X), which will be clearly visible in the redshift dependence of dN /dX, the 
zeroth moment of f(N,X), and of £l g (z), the first moment of f(N,X) (see § 
2.3 and § 2.4). At low column densities, f(N,X) increases with redshift by a 
factor of 2 at z > 2.2. Prochaska, Herbert-Fort, & Wolfe (2005) detect a similar 
evolution at higher values of N. By contrast the shape of f(N,X) does not 
appear to evolve with redshift. This is in disagreement with the earlier results 
of Storrie-Lombardi & Wolfe (2000) and Peroux et al. (2003b) who used much 
smaller samples to claim that f(N,X) steepened at z > 3.5.. 

2.3 dU/dX 

In Figure 4 we plot the most recent evaluation of dN / dX versus z for damped 
Lya systems with z>0 (see Prochaska, Herbert-Fort, & Wolfe 2005). The solid 
line traces the value of dM jdX derived in a series of 0.5 Gyr time intervals to 
reveal the effects of binning. The data points at z = are three estimates of 
{dN /dX) z= o based on the H I properties of nearby galaxies. The figure shows a 
decrease by a factor of two in dN jdX from z=4 to 2. From equations 2 and 3 
we see that the decrease in dN jdX reflects a decrease in either H I cross section, 
comoving density, or of both quantities. Prochaska, Herbert-Fort, & Wolfe (2005) 
use the Press-Schechter formalism to show that significant variations in comoving 
density with time are unlikely to occur. Therefore, within the context of CDM 
models the most likely explanation for the changes in dN /dX is a decrease in H 
I cross section with time. This is probably due to feedback mechanisms such as 
galactic winds. 

Figure 4 also shows that dN / dX at z « 2 is consistent with the present-day 
value; i.e., the data are consistent with an unevolving population of galaxies 2 . 
By comparison, Wolfe et al. (1986) required more than a factor of 4 increase with 
redshift in dN jdX. The discrepancy arises from the earlier use of an Einstein- 
deSitter rather than the modern ACDM cosmology to estimate AX intervals, 
and from the lower values of (dN /dX) z= o used in the earlier work. This result 
implies that between z = 1 and 2 smaller galaxies merged to produce bigger ones 
such that the product of comoving density and total H I cross-section for iV(H I) 
> 2xl0 20 cm -2 is conserved. 

2.4 n g (z) 

Figure 5 shows the most recent determination of £l g (z). From Equation (4) we 
see that £l g (z) is sensitive to the upper limit N max unless a « — 2. This led 
to large uncertainties in £l g (z) in previous work because a was not measured 

2 Ryan- Weber, Webster, & Staveley-Smith (2003) present evidence that f(N, X) at z—Q is 
significantly lower in amplitude than the results at higher redshifts, but this result is puzzling 
since comparison between the resultant dAf/dX at z—Q with the higher redshift data in Figure 4 
reveals no evidence for evolution. 
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with sufficient accuracy to rule out a > — 2. However, with the large sample 
of over 600 damped Lya systems Prochaska, Herbert-Fort, Sz Wolfe (2005) use 
several tests to show that £l g (z) converges. First, they compute a\ for a single 
power-law fit to f(N,X) by increasing N m i n from 2xl0 20 cm -2 . Using the full 
sample of damped Lya systems they find a,\ decreases with increasing N m i n from 
-2.2 at iV m j n =2xl0 20 cm~ 2 to less than -3 at N rnin > 10 21 cm -2 . At the same 
time they find that ct\ is insensitive to variations in N max . Second, they compute 
the sensitivity of £l g (z) to N max . Both the double power-law and T function 
solutions converge to the value indicated by the data (Equation 5). By contrast 
the single power-law solution does not converge. This is the first evidence that 
£lg(z) converges by iV ~ 10 22 cm -2 . 

Next consider the redshift evolution of Q g (z). Starting at the highest redshifts, 
no increase of £l g {z) with decreasing z is present at z > 3.5, contrary to earlier 
claims (Storrie-Lombardi & Wolfe 2000; Peroux et al. 2003b). On the other 
hand, Figure 5 shows the first statistically significant evidence that £l g (z) evolves 
with redshift. Specifically, Q, g {z) decreases from lxlO -3 at z=3.5 to 0.5 xl0~ 3 
at z=2.3, which mirrors the decline in dM jdX discussed in § 2.3. The same 
mechanism is likely to cause the decline in both quantities; i.e., a decrease in H I 
cross section due to feedback. But at z < 2.3 the picture is somewhat confusing. 
Figure 5 shows an increase of £l g (z) by z ~ 2, which is consistent with the values 
of £l g (z) in the two lower redshift bins at < z < 2. Indeed, the data are 
consistent with no evolution, if one ignores the redshift interval centered at z 
= 2.3. However, Prochaska, Herbert-Fort, and Wolfe (2005) emphasize that the 
uncertainties in the data at < z < 2.3 are much larger than at z > 2.3, and 
thus such conclusions should be treated with caution. 

Next, we compare the high-z values of Q g (z) with various mass densities at 
z = 0. First, comparison with the current density of visible stars, Q*(z = 0), 
reveals that Q g {z) at z ~ 3.5 is a factor of 2 to 3 lower than Q*(z = 0): if the 
census of visible stars were restricted to stellar disks, then fl g (z) at these redshifts 
would exceed = 0). A straightforward interpretation of this concurrence is 
that damped Lya systems provide the neutral gas reservoirs for star formation at 
high redshifts. However, since J7*(z = 0) exceeds f2 fl (z~3.5), the reservoir must 
be replenished with new neutral gas before the present epoch. Further evidence 
for replenishment is that it is required to compensate for gas depletion due to 
star formation detected in damped Lya systems (see § 8). As a result the "closed 
box" hypothesis for evolution in damped Lya systems is unlikely to be correct 
(see Lanzetta et al. 1995). 

Figure 5 also shows that fl g (z) at z ~ 3.5 is significantly higher than £l g (z) 
at z = 0, which is deduced from surveys for 21 cm emission. Therefore, damped 
Lya systems provide direct evidence for the widely held theoretical view that the 
neutral gas content of the Universe was larger at high redshifts than it is today. 
Figure 5 also shows that £l g (z) at z ~ 3.5. is at least a factor of 10 greater than 
= 0) in dwarf irregular galaxies, which argues against the idea that damped 
Lya systems evolve into such objects (e.g. Jimenez et al. 1999). 

Therefore, since Lyman limit systems do not contribute significantly to the 
neutral gas content at any redshift (see Prochaska, Herbert-Fort, & Wolfe 2005) 
and ignoring the possible existence of a significant population of dusty giant 
molecular clouds, we conclude that damped Lya systems contain most of the gas 
available for star formation at z > 1.6. At low redshifts the ionizing background 
is reduced and lower N(K I) systems might be mainly neutral. But at z = 0, 
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Minchin et al. (2003) find a paucity of galaxies with column densities less than 
iV(H I) = 2x 10 20 cm -2 measured from 21cm emission, implying that at the lowest 
redshifts, damped Lya system column densities comprise most of the neutral gas 
in the Universe. As a result, damped Lya systems dominate the neutral-gas 
content of the Universe in the redshift interval z = [0,5]. 

Of course, all of these conclusions ignore obscuration by dust in damped Lya 
systems, which may have biased the form of f(N,X) (see Pei, Fall, & Bechtold 
1991; Fall & Pei 1993). We discuss this possibility in § 10. They also ignore 
biasing due to lensing, which may be present (see § 11). 

2.5 Model Comparisons 

Here we discuss attempts to use the H I content of damped Lya systems to test 
models of galaxy formation and evolution. 

2.5.1 Passive Evolution: The Null Hypothesis 

The evolution of H I content within the null hypothesis of passive evolution 
has been modeled by Boissier, Peroux, & Pettini (2003). In this scenario damped 
Lya absorption arises in disk galaxies with the comoving density of current spiral 
galaxies. The models are hybrids of passive evolution, in which the H I content is 
changed only by processes of stellar evolution, and the spherical collapse model in 
an expanding Universe, in which high-z disks are smaller than current disks. The 
models are successful in explaining the lack of evolution in dAf/dX, f(N, X), and 
Qg(z) at z < 2 but fall short at higher redshifts because of delayed disk formation. 
For these reasons, the authors suggest the added presence of a population of low- 
surface brightness, gas-rich galaxies at z > 2. However, evidence against the 
"closed box" hypothesis discussed in § 2.4 is a difficult challenge for this and all 
passive evolution models. 

2.5.2 Numerical Simulations 

Several aspects of cosmology and galaxy formation have been examined through 
comparisons of numerical simulations of galaxy formation in a CDM Universe 
with the observed HI properties of damped Lya systems. A critical feature of 
these models is that gas falling onto dark-matter halos is heated to their virial 
temperatures, then cools off, and collapses toward the central regions of the halos. 
Galaxies arise from the formation of stars out of the cool (presumably neutral) 
collapsed gas and evolve through mergers between dark-matter halos and further 
infall of gas onto the halos. 

The first studies (Ma & Bertschinger 1994; Klypin et al. 1995) constrained the 
cosmological mass constituents through comparisons with £l g (z). The observa- 
tions severely restricted the contribution from a hot component (i.e. neutrinos) 
as these cold+hot dark matter cosmogonies underpredicted structure formation 
at early times (Katz et al. 1996). Subsequent papers by Gardner et al. (1997a,b, 
2001) examined the properties of damped Lya systems in their smooth parti- 
cle hydrodynamic (SPH) simulations. These models generically underpredicted 
the incidence of damped Lya systems, which the authors argued was due to 
an insufficient mass resolution of W u Mq. They did find reasonable agreement 
with the data, however, by extrapolating to halos with > 10 10 M using the 
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Press-Schechter formalism, and by assuming that the HI cross-section followed 
the power-law expression A oc v x c % . As stressed by Prochaska & Wolfe (2001), 
this power-law expression implies a v c distribution that is incompatible with the 
observed damped Lya velocity widths (see § 6). 

Nagamine, Springel, & Hernquist (2004a) recently analyzed a comprehensive 
set of high resolution SPH simulations of a ACDM Universe. By contrast with 
Gardner et al. (1997a) they find that halo masses down to Mh ~ 10 8 M@ con- 
tribute to the H I cross-sections: halos with Mh < 10 8 M Q do not contribute 
since they contain only photo-ionized gas. In turn, Nagamine et al. find a 
steeper power-law expression A oc v^ 7 , i.e., massive halos make a larger contri- 
bution to dhf jdX than in previous models. They conclude that Gardner et al. 
(1997a) predicted an overabundance of damped Lya systems with < 10 10 
M Q because the slope of their A versus v c relation was too shallow. Nagamine, 
Springel, & Hernquist (2004a) were the first to include mass loss of neutral gas 
due to winds, which increases the median halo-mass contribution to dM jdX to 
10 11 Mq. Winds also prevent an overabundance of Q g (z) at z > 4, but under- 
predict Q g (z) at z < 4. Nagamine et al. (2004a) continue to find a deficit of 
damped Lya systems with iV(H I) < 10 21 cm -2 . The deficit of systems with low 
7V(H I) is a generic effect seen in most (e.g. Figure 3 in Katz et al. 1996) but not 
all numerical simulations (Cen et al. 2003) and is a shortcoming that needs to 
be addressed. On the other hand, the Nagamine, Springel, & Hernquist (2004a) 
models are the most successful in reproducing the evolution of tt g (z)) at z > 2 
(see Figure 5). 

2.5.3 Semi-analytic and Analytic Models 

The semi-analytic models were proposed to include processes beneath the res- 
olution of the numerical simulations with phenomenological descriptions of star 
formation, gas cooling, and the spatial distribution of the gas. The latter is 
included since the simulations failed to reproduce the correct sizes and angu- 
lar momenta of present-day galactic disks (Navarro & Steinmetz 2000), and we 
do not know whether the simulations produce the correct spatial distribution of 
neutral gas at z ~ 3. This is a concern because damped Lya systems are a cross- 
section weighted population of high-redshift layers of neutral gas and therefore 
the results will be sensitive to the gas distribution at large impact parameters. 
By contrast with the numerical simulations, one uses analytic expressions from 
Press-Schechter theory or its extensions (Press & Schechter 1974;Sheth, Mo, & 
Tormen 2001) to compute the mass function of halos that evolves from a given 
power spectrum. 

Mo & Miralda-Escude (1994) modeled damped Lya systems with the Press- 
Schechter formalism in both mixed Cold+Hot dark matter and ACDM cosmolo- 
gies. Kauffmann (1996) used a SCDM cosmology ([n M ,^A,^]=l-0,0.0,0.5) to 
construct improved semi-analytic models for damped Lya systems. She assumed 
spherical geometry for the halos and let the neutral gas in a given halo be confined 
to a smaller, centrifugally supported disk. To compute the radial distribution of 
the neutral gas, the angular momentum per unit mass of the disk was set equal 
to that of the halo. Using Monte Carlo methods, she computed the formation 
and growth of individual halos with time. The neutral gas content was assumed 
to be regulated by accretion due to mergers and star formation, but feedback due 
to winds was omitted. In common with the Nagamine, Springel, & Hernquist 
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(2004a) simulations, the Kauffmann (1996) models (1) reproduced the Q g (z) re- 
lation, and (2) exhibited a deficit of systems with iV(H I) < 10 21 cm -2 . 

Mo, Mao, & White (1998) constructed models for disk formation that were 
also based on the Press-Schechter formalism. These models extended the work 
of Kauffmann (1996) by considering disks drawn from a distribution of halo spin 
parameters, Xh, rather than using Kauffmann's technique of assigning the mean 
value of Xh to each disk. Since disks detected in a survey for damped Lya sys- 
tems are drawn from a cross-section weighted sample favoring bigger disks, the 
distribution of Xh will be skewed to values higher than the unweighted mean. 
The result is larger H I cross-sections and higher detected rotation speeds. Con- 
sequently, they found agreement between the predicted and observed dhf jdX 
relation. Mailer et al. (2001) then suggested a model in which the gas is in ex- 
tended Mestel (1963) disks in which the surface density falls off inversely with 
radius. In this case the disks overlap and as a result the observed f(N,X) is re- 
produced; i.e., there is no deficit of systems with iV(H I) < 10 21 cm~ 2 . However, 
it is unclear whether such extended disks will either form or survive sufficiently 
long to contribute to the H I cross-sections of damped Lya systems. Furthermore, 
most semianalytic models overestimate Q g (z) at z > 2, since they underpredict 
feedback processes at these redshifts (see Figure 5). 

3 CHEMICAL ABUNDANCES 

Because the damped Lya systems comprise the neutral gas reservoir for star 
formation at high redshifts, a determination of their metal content is a crucial 
step for understanding the chemical evolution of galaxies. Therefore, the mass 
of metals per unit comoving volume that they contribute indicates the level to 
which the neutral gas reservoir has been chemically enriched. Since the metal 
abundances of damped Lya systems have been determined in the redshift interval 
z = [0, 5], it is now possible to track the chemical evolution of the reservoir back 
~ 10 Gyr to the time the thin disk of the Galaxy formed [z = 1.8 for the WMAP 
cosmological parameters adopted here), and to earlier epochs. As a result, one 
can construct an "age-metallicity" relation not just for the solar neighborhood 
(see Edvardsson et al. 1993) but for a fair sample of galaxies in the Universe. 

In this section we describe the main results that have emerged from abundance 
studies of damped Lya systems. This subject has recently been reviewed in an 
excellent article by Pettini (2004). 

3. 1 Methodology 

The element abundances of the damped Lya systems are the most accurate mea- 
surements of chemical enrichment of gas in the high-redshift Universe. The mea- 
surements are accurate for several reasons: (1) For the majority of damped Lya 
systems, hydrogen is mostly neutral, i.e., H°/H=l, and most of the abundant el- 
ements are singly ionized, though a minority are neutral, i.e., Fe + /Fe=Si + /Si=l, 
etc., while O /O=N°/N=l. The singly ionized elements have ionization poten- 
tials of their neutral states that are lower than the ionization potential of hydro- 
gen, IP(H) (=13.6 eV). With Lyman limit optical depths, tll » 10 3 , damped 
Lya systems are optically thick at photon energies, IP(H) < hv < 400 eV. As a 
result, only photons with hv < IP(H) and hv > 400 eV penetrate deep into the 
neutral gas. When FUV photons with hv < IP(H) penetrate, they photoionize 



14 



Wolfe, Gawiser & Prochaska 



the neutral state of each element to the singly ionized state. But this state is 
shielded from photons with IP(H) < hv < 400 eV, which would otherwise pho- 
toionize the elements to higher states. Photons with hv > 400 eV will produce 
species that are doubly ionized (e.g. Fe ++ and Al ++ ) and singly ionized (e.g. 
Ar + ), but because the photoionization cross-sections are low at such high photon 
energies, the ionization rates are low. It is possible to detect all of these species 
because they exhibit resonance transitions that are redshifted to optical wave- 
lengths accessible with ground based spectrographs. (2) In § 1 we saw that Voigt 
fits to the damped Lya profiles result in typical errors of 0.1 dex in iV(H I). As we 
shall see, the errors in the column densities which give rise to the narrow low-ion 
lines are typically 0.05 dex. Consequently, errors in [X/H] 3 are relatively low 
(typically about 0.1 dex). (3) Column densities are straightforward to measure 
from resonance lines, since their optical depths are independent of poorly deter- 
mined physical parameters such as the density and temperature of the absorbing 
gas. 

By contrast, abundance determinations for the other constituents of the high- 
redshift Universe are more uncertain primarily because the gas is ionized. As 
a result, the abundances are subject to ionization corrections which depend on 
uncertainties in the shape of the ionizing continuum radiation and on the trans- 
port of such radiation. Furthermore, the strengths of QSO emission lines depend 
on the temperature and density of the emitting gas as well as uncertain photon 
escape probabilities in the case of resonance scattering. Typical error estimates 
are about 50 % per object, which is several times higher than for damped Lya 
systems (see Hamann & Ferland 1999 for an excellent review of this subject). 

Figure 6 shows examples of absorption profiles obtainable with the HIRES 
Echelle spectrograph mounted on the Keck I 10-m telescope. The figure shows 
velocity profiles for abundant low ions in two damped Lya systems. As in most 
damped Lya systems the gas that gives rise to low-ion absorption lines in these 
two objects is comprised of multiple discrete velocity structures of enhanced den- 
sity; i.e., clouds. To infer the ionic column densities required for element abun- 
dance determinations, one integrates the "apparent optical depth" (Savage & 
Sembach 1991; Jenkins 1996) over the velocity profile. 

To illustrate the essentials of abundance determinations we focus on two damped 
Lya systems, one metal-poor (DLA1108— 07 at z = 3.608) and the other metal- 
rich (DLA0812+32 at z = 2.626) 4 . The corresponding velocity profiles in Fig- 
ure 6 describe the challenges as well as the advantages of measuring damped 
Lya abundances. First consider the challenges. The abundance of carbon has 
not been accurately determined for any damped Lya system because the only 
resonance transition outside the Lya forest, C II A 1334.5, is not only saturated 
in the metal-rich system (Figure 6b), but is also saturated in the metal-poor 
system (Figure 6a). The availability of several O I transitions makes it possible 
to place bounds on the oxygen abundance. Several authors used saturated O 
I A 1302.1 to obtain lower limits and the weaker O I A 971.1 or O I A 950.8 
transitions for upper limits since the latter transitions are usually blended with 
Lya forest absorption features (D'Odorico & Molaro 2004; Molaro et al. 2000; 
Dessauges-Zavadsky, Prochaska, k, D'Odorico 2002). On the other hand, a direct 

3 Here and in what follows the relative abundance of elements X and Y is defined with respect 
to the solar abundance on a logarithmic scale; i.e. [X/Y]= logio(X/Y)— logio(X/Y)0. 

4 Here and in what follows we designate a damped Lya system toward a QSO with coordinates 
hhmmideg as DLAhhmmideg 
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determination of [0/H] is possible for the metal-rich system shown in Figure 6b 
because this is the only known case in which an unsaturated transition, O I] A 
1355.6, is detected. 

Second consider the advantages. Abundance determinations are possible for 
Fe, Si, and S because of the presence of transitions with a wide range of oscillator 
strengths. In the case of Fe, the oscillator strengths /1611.2 = 0.00136 for Fe II A 
1611.2 and /i 6 os.4 = 0.0580 for Fe II A 1608.4. Thus, in the metal-poor damped 
Lya system in Figure 6a Fe II A 1611.2 is undetected, while unsaturated Fe II 
A 1608.4 is detected. By contrast, in the metal-rich damped Lya system in Fig- 
ure 6b unsaturated Fe II A 1611.2 is detected, while Fe II A 1608.4 is saturated. 
In both systems the iron abundance is determined from the unsaturated transi- 
tions. Similarly, Figure 6 also demonstrates how the Si II A 1304.3, 1808.0 pair of 
transitions determines the silicon abundance for the two damped Lya systems. 



3.2 Metallicity 

Whereas the abundance ratios discussed above refer only to elements in the gas 
phase, some fraction of each element could be depleted onto dust grains, as in the 
Galaxy ISM (Jenkins 1987). This possibility was recognized early on by Meyer, 
Welty, & York (1989) and Pettini, Boksenberg, & Hunstead (1990) who made use 
of the Zn II AA 2026.1, 2062.6 doublet to measure the metallicities of damped Lya 
systems. Zn is well suited for this purpose because Zn is relatively undepleted 
in the ISM with a mean depletion of [Zn/H] pa —0.23 (Savage & Sembach 1996). 
Moreover Zn was believed to be an accurate tracer of Fe peak elements since 
[Zn/Fe] pa for stars with metallicities, —2.0 < [Fe/H] < (but see discussion 
below). In addition the combination of the low solar abundance of Zn and the 
oscillator strengths of the Zn II transitions implies they should be unsaturated 
for iV(H I) < 10 21 cm -2 , provided the velocity dispersion of the g > 4 

km s -1 . Because of its proximity in wavelength, the Cr II AA 2056.2, 2062.2, 
2066.1 triplet was used to study depletion, since most of the Cr in the ISM is 
locked up in grains (Jenkins 1987). 

In subsequent surveys on several 4 m class telescopes, Pettini and colleagues 
(Pettini et al. 1994, 1997b, 1999) increased the size of their sample and confirmed 
that damped Lya systems are metal-poor in the redshift interval z = [0.5, 3.0]. 
Pettini (2004) found that the cosmic metallicity < Z > = — 1.11±0.38, where 
< Z > is defined as the log of the ratio of the comoving densities of metals and 
gas, Q m etals/Qg, relative to the solar abundance; i.e. from Equation 5 



< Z >= log 10 



£"=1 lOW^JVj 



-log lo (M/H) , (6) 



where M stands for the metallicity indicator, which in this case is Zn. Second, 
surprisingly, there is no positive evidence for redshift evolution. Specifically, 
Pettini (2004) finds no statistically significant evidence for redshift evolution in 
< Z >. This is contrary to most models of chemical evolution (see § 9), which 
predict an increase in the mean Zn abundance with decreasing redshift, and 
further predict that the metallicity should approach < Z >= 0, by the current 
epoch. The sub-solar values of < Z > at z < 1 raised the possibility that damped 
Lya systems do not evolve into normal current galaxies (Pettini et al. 1999). 

Further progress was achieved with the completion of a larger survey of over 
120 damped Lya systems carried out primarily on the Keck 10-m telescopes 
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(Prochaska et al. 2003a). In this survey most of the metallicities, [M/H], are 
obtained from measurements of a-enhanced elements Si, S, and O in order of 
decreasing priority and in a few cases from Zn. Like Zn, S and O are volatile 
elements which are essentially undepleted in the ISM. While the refractory ele- 
ment Si is depleted in the ISM, it is only mildly depleted in damped Lya systems, 
where Si tracks S, i.e., [Si/S] > -0.1 (Prochaska & Wolfe 2002), and thus can 
generally be used as an unbiased metallicity tracer. Furthermore, since S and 
Si have higher solar abundances than Zn, they can be used to probe down to 
metallicities below the Zn threshold of [Zn/H]« —1.7. In addition, the shorter 
wavelengths of crucial transitions such as S II A 1250.5 and Si II A 1304.3 allow 
one to obtain metal abundances at higher redshifts than are accessible with the 
Zn II transitions alone. Note, the idea of combining abundances of Zn and a- 
enhanced elements is plausible if Zn is a tracer of elements such as S, Si, etc. The 
recent finding by Prochaska & Wolfe (2002) that [Si/Zn] = 0.03±0.05 supports 
this hypothesis. Further support comes from the finding that [Zn/Fe] ranges 
between 0.10 to 0.20 (Prochaska et al. 2000; Nissen et al. 2004; Bihain et al. 
2004) in stars with [Zn/H] > —1.5, which indicates Zn is not a strict tracer of Fe 
peak elements. In fact there is currently little reason to expect [Zn/Fe] =0 aside 
from a coincidence related to the star formation history of the Galaxy (Fenner, 
Prochaska, & Gibson 2004). 

The results of Prochaska et al. (2003a) are shown in Figure 7 (updated to 
include new data at z < 1.5 from Kulkarni et al. 2005 and Rao et al. 2005). 
The new survey confirms the low metallicities of damped Lya systems found by 
Pettini and colleagues. However, the greater accuracy and larger redshift range 
of the new survey allows one to draw additional conclusions. First, there are no 
damped Lya systems with [M/H] < —2.6. This limit is robust because there are 
no damped Lya systems without significant metal absorption. Second, Prochaska 
et al. (2003a) find statistically significant evidence for a linear increase of < Z > 
with decreasing z. This result is robust owing to the large value of J2T ^i- This 
is important since the shape of f(N,X) indicates that < Z > is sensitive to the 
metallicity of systems with the largest values of iV(H I) . Because J2T > 1 x 10 22 
cm~ 2 in each of the high-redshift bins, only unusual, very metal-rich systems with 
7V(H I) > 10 22 cm~ 2 could increase < Z > significantly; i.e., only systems which 
depart significantly from the current iV(H I) versus [M/H] relation could cause a 
marked increase in < Z > . Earlier claims for evolution had statistical significance 
lower than 3a and sampled lower values of J2? Ni (Kulkarni & Fall 2002; Vladilo 
2002b). 

The "age-metallicity" relationship depicted in Figure 7 provides new infor- 
mation about the enrichment history of damped Lya systems. Specifically, the 
absence of any system with a metallicity [M/H] < — 2.6 sets the damped Lya 
systems apart from the Lya forest. From their analysis of the Lya forest, Simcoe, 
Sargent, & Rauch (2004) find a median abundance, [C,0/H] = —2.8 and find that 
30 % of their sample have [C,0/H] < —3.5. Schaye et al. (2003) find similar re- 
sults for [C/H]. While they deduce a higher median abundance for Si, i.e., [Si/H] 
= —2.0, about 40 % of their systems are predicted to have [Si/H] <— 2.6. Clearly 
the bulk of the damped Lya population has a different enrichment history than 
the Lya forest. To explain the presence of the metallicity floor, Qian & Wasser- 
burg (2003) use a standard chemical evolution model to show that star formation 
in damped Lya systems results in a rise in metal abundance which is so rapid 
that the probability for detecting systems with [M/H] < — 2.6 is exceedingly 
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small. 

Figure 7 also poses several dilemmas for models of chemical evolution. First, 
if most of the gas in damped Lya systems in the redshift interval z = [1.6, 
4.5] were converted into stars, then most of the stellar mass in current galaxies 
would be metal poor, contrary to observations Tremonti et al. (2004). Second, 
the age-metallicity relation of the thin disk of the Galaxy (Edvardsson et al. 
1993), indicates that the thin disk formed at lookback times less than 10 Gyr 
(i.e., z f» 1.8) and that chemical enrichment proceeded such that all thin disk 
stars formed with [M/H] > —1.0. But the lower panel in Figure 7 shows that 
[M/H] < —1.0 in about half of the damped Lya systems with look-back times 
under 10 Gyr. While this result is subject to the uncertainties of small number 
statistics and observational bias, the current metallicity trends in low-redshift 
damped Lya systems suggest that damped Lya systems may not trace the star 
formation history of normal galaxies (Pettini et al. 1999). Third, if the linear 
increase of < Z > with decreasing redshift deduced at z > 1.6 is extrapolated to 
z = 0, the current mean metal abundance of galaxies would be equal to —0.69 
which appears too low. But since the age-metallicity relationship is essentially 
unconstrained by the data at z < 1.6, such extrapolations should be treated with 
caution. Indeed, < Z > is doubling every Gyr at z > 2, and if we assumed 
< Z > to be a linear function of time rather than redshift, then we would find 
that < Z >«0 by z « 0.5. 

3.3 Relative Abundances 

In § 3.2 we described evidence that damped Lya systems are metal-poor. We 
discussed measurements of Zn and Cr which indicate a gas-phase abundance 
ratio, [Zn/Cr] > 0, implying depletion of Cr by dust. Since metal-poor stars 
in the Galaxy exhibit different nucleosynthetic abundance ratios than the sun 
(Wheeler, Sneden, & Truran 1989), the abundance patterns observed in damped 
Lya systems are probably due to some combination of nucleosynthetic and dust 
depletion patterns. In this section we briefly describe efforts to unravel these 
effects. The reader is referred to a series of recent papers for a more thorough 
discussion of these issues (Prochaska & Wolfe 2002; Vladilo 2002b; Pettini 2004). 

3.3.1 Depletion 

The discussions of metallicity in the previous subsections implicitly assumed that 
deviations of (X/H) ffas , the gas- phase abundance of element X, from the so- 
lar abundance, (X/H)©, were only due to changes in the intrinsic abundance, 
(X/H)j nt . However, as mentioned previously, (X/H) gas will also deviate from 
(X/H) Q if element X is depleted onto grains. One of the major challenges in 
damped Lya research is to untangle these two effects. 

The traditional method used by most workers in the field is to compare the 
abundance of refractory element, X, to volatile element, Y, for which (X/Y)j n j 
= (X/Y) in stars with a wide range of absolute abundances. In that case the 
condition [X/Y] 7^ is unlikely to have a nucleosynthetic origin. Rather it likely 
arises from depletion of the refactory element onto grains. Such a comparison is 
made in Figure 8a, which is a plot of [Zn/Fe] versus [Zn/H] for a sample of 32 
damped Lya systems. The figure reveals an unambiguous correlation between 
[Zn/Fe] and [Zn/H]: a Kendall r test rules out the null hypothesis of no cor- 
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relation at more than 99.7 % confidence. Because [Zn/Fe] < 0.2 for Galactic 
stars with [Fe/H] > —2.0 , the most plausible explanation for this correlation is 
that in damped Lya systems the depletion of Fe onto grains increases with metal 
abundance. This argument also suggests that the depletion level decreases with 
decreasing metal abundance. In that case [Zn/Fe] should approach the intrinsic 
nucleosynthetic ratio, [Zn/Fe]j n t, in the limit [Zn/H] << 0. Determination of 
[Zn/Fe] j n £ is important as it indicates the nucleosynthetic history of these ele- 
ments (Hoffman et al. 1996), and it is required for determining the dust-to-gas 
ratio, k. For example, Wolfe, Prochaska, h Gawiser (2003b) show that 



where in this case X=Fe and Y=Zn. 

Our discussion emphasizes the importance of estimating the intrinsic, nucle- 
osynthetic ratio, [Zn/Fe]j„(. On the other hand, the observed Zn abundances 
are not sufficiently low for the asymptotic approach to [Zn/Fe] i n t to be detected. 
Specifically, because the Zn II transitions are weak, only two damped Lya sys- 
tems with [Zn/H] < —1.5 have been detected (Lu, Sargent, & Barlow 1998; 
Molaro et al. 2000; Prochaska & Wolfe 2001). By contrast, clouds of such low 
metallicity can be easily detected in the strong Si II transitions, as shown in Fig- 
ure 8b, which plots [Si/Fe] versus [Si/H] down to [Si/H] = —2.6. The figure gives 
convincing evidence that in the limit of vanishing metallicity, [Si/Fe] approaches 
~ 0.3 rather than 0. Furthermore, at metallicities [Si/H] > —1 we see evidence 
for an increase in [Si/Fe] with increasing [Si/H]. This is the same phenomenon 
seen in the [Zn/Fe] versus [Zn/H] diagram, which we plausibly attributed to dust. 
The amplitude of the increase is weaker for [Si/Fe] because Si is weakly depleted. 
On the other hand the increase of [Si/Fe] with [Si/H] is stronger evidence for 
dust since the nucleosynthetic origin of Si is better understood than that of Zn 
(Hoffman et al. 1996)). 

3.3.2 Nucleosynthetic Abundance Patterns 

• a Enhancements? 

In § 3.3.1 we argued that the asymptotic behavior exhibited by the [Si/Fe] ratio 
in the limit [Si/H] << (Figure 8b) indicated a nucleosynthetic ratio, [Si/Fe] j n t 
~ 0.3. This asymptotic limit is robust, as it is based on a large number (56) 
of precision measurements obtained with echelle spectrometers on 8- to 10-m 
class telescopes. It also has important implications for the chemical evolution of 
damped Lya systems if it equals the intrinsic nucleosynthetic ratio. The reason 
is that disk stars in the Galaxy exhibit a systematic decrease of [a/Fe] with 
increasing [Fe /H] , which indicates the increase with time of Fe contributed to the 
Galaxy ISM by type la supernovae relative to type II supernovae (Edvardsson 
et al. 1993). The presence of such trends in damped Lya systems would support 
the argument that they are the progenitors of ordinary galaxies. However, the 
existence of intrinsic a enhancements in damped Lya systems is controversial. 
Using the Vladilo (1998, 2002a) models, Vladilo (2002b) and Ledoux, Bergeron, 
&; Petitjean (2002) examined [Si/Fe] ratios corrected for depletion effects and 
found median values of [Si/Fe] j n t consistent with solar. Similarly, several studies 
of the depletion-free [O/Zn] and [S/Zn] ratios resulted in [a/Zn] ratios below 
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those of metal-poor stars (Molaro et al. 2000; Centurion et al. 2000; Nissen et al. 
2004). 

Is it possible to resolve these conflicts? The lower values of [a/Zn] are compati- 
ble with the higher value of [a/Fe]^ indicated by Figure 8b since [a/Fe] = [a/Zn] 
+ [Zn/Fe] and Prochaska et al. (2000) and Chen, Kennicutt, & Rauch (2005) 
find [Zn/Fe]«0.15 for thick disk stars with -0.9 < [Fe/H] < -0.6, while Nissen 
et al. (2004) find [Zn/Fe] « 0.1 for stars with [Fe/H] <-1.8. Both results are 
consistent with [Si/Fe]j„i~ 0.2 to 0.4. If such a enhancements are confirmed in 
damped Lya systems, one would conclude that the depletion corrections used by 
Vladilo (2002a) and Ledoux, Bergeron, & Petitjean (2002) were too large. The 
latter are compatible with the dust content suggested by the Pei, Fall, & Bechtold 
(1991) study of reddening in damped Lya systems. But, since the more recent 
study of Murphy & Liske (2004) argues against such a high dust content, the 
depletions used to correct the [Si/Fe] ratios may be too large. We also note that 
[a/Fe}i n t=0 would imply significant depletion of Fe at [Si/H] < —1, which is 
apparently at odds with the insensitivity of [Si/Fe] to increases in [Si/H] shown 
in Figure 8b. But this behavior may result from two compensating effects: an 
increase in [Si/Fe] due to Fe depletion and a decrease in [Si/Fe] due to increasing 
Fe enrichment from type la supernovae. Because of these uncertainties, it may 
be premature to use the [a/Fe] ratios in damped Lya systems as discriminants 
between competing galaxy formation scenarios (e.g. Tolstoy et al. 2003; Venn 
et al. 2004). 

• Nitrogen Enrichment 

Pettini, Lipman, & Hunstead (1995) first detected nitrogen in damped Lya 
systems and suggested the [N/a] versus [a/H] plane could be used as a clock to 
infer their ages. According to Henry, Edmunds, & Koppen (2000) a burst of star 
formation would coincide with the injection of a elements into the surrounding 
ISM by type II supernovae, followed by the injection of 14 N by AGB stars more 
than 0.25 Gyr later. In the local Universe, one identifies a plateau of [N/a] values 
(with value ps — 0.7dex) at low metallicity presumably consisting of objects with 
ages greater than 0.25 Gyr. Within this interpretation, metal-poor objects with 
ages less than 0.25 Gyr would have systematically lower [N/a] values while more 
evolved damped Lya systems, in which 14 N production has caught up, would have 
[N/a]«-0.7. Recent studies (Prochaska et al. 2002a; Pettini et al. 2002; Cen- 
turion et al. 2003) have shown that the majority of [N/a] values for the damped 
Lya systems are near the plateau but there is a population of damped Lya sys- 
tems with [N/a] —1.5 and very few damped Lya systems with intermediate 
[N/a] values. The damped Lya systems with [N/a]=— 0.7 must be older than 
0.25 Gyr, indicating they are not transient objects as suggested in some schemes 
(see Qian & Wasserburg 2003) but rather have ages comparable to the age of the 
Universe at z ~ 3. 

The observations also pose a challenge to interpreting the [N/a] value as a 
strict age diagnostic. If the ages of the damped Lya systems are comparable 
to 2.5 Gyr, the age of the Universe at z ~ 2.5, then fewer than 10% of the 
objects would have [N/a]=— 1.5, contrary to current observation. Prochaska et al. 
(2002a) interpret the paucity of systems with —1.5 < [N/a] < —0.7 as evidence 
for a bimodal IMF where systems near the plateau at [N/a]=— 1.5 are drawn 
from an IMF truncated from below at M=7.5 M Q . In this scenario, damped 
Lya systems near the plateau at [N/a]=— 1.5 need not be younger than 0.25 
Gyr, while systems near the plateau at [N/a]=— 0.7 are objects older than 0.25 
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Gyr in which N production is due to the full mass range of intermediate-mass 
stars drawn from a standard IMF. More recently Molaro (2003) argued against 
a bimodal IMF by suggesting that damped Lya systems near the [N/a]=— 1.5 
plateau are younger than the 0.25 Gyr "catch-up" time. Meynet & Maeder (2002) 
and Chiappini, Matteucci, & Meynet (2003) suggested a mechanism for obtaining 
a more uniform distribution between the two plateaus. They showed that stellar 
rotation causes enhanced mixing between the H-burning and He-burning layers, 
thereby producing greatly enhanced 14 N production in massive stars. Meynet 
& Maeder (2002) reproduced the [N/a]=-1.5 plateau for stars withM=8-120 
Mq for rotation speeds fsim=400 km s _1 . Moreover, rotation may extend the 
effective lag time between N and a production for intermediate mass stars beyond 
0.25 Gyr. However, the lack of many damped Lya systems with [N/a] in between 
— 1.5 and — 0.7dex argues against this mechanism and in favor of the bimodal 
IMF. 

• Metal-Strong damped Lya system 

There exists a small subset of damped Lya systems for which the product 
of HI column density and metallicity imply very strong metal-line transitions. 
These 'metal-strong' damped Lya systems yield abundance measurements for 
over 20 elements including O, B, Ge, Cu, and Sn. Figure 9 shows the elemental 
abundance pattern obtained for DLA0812+32, the z = 2.626, metal-rich damped 
Lya system discussed in § 3.1 (Prochaska, Howk, & Wolfe 2003c). This is the 
first damped Lya system for which absolute abundances for B and Ge have been 
determined, and it is one of the few objects for which an accurate measurement 
of [O/H] is possible. With the detection of over 20 elements, the metal-strong 
damped Lya systems permit a global inspection of its enrichment history. The 
dotted line in Figure 9 is the solar abundance pattern scaled to oxygen. The good 
match to the data shows that this system exhibits an enrichment pattern resem- 
bling that of the Sun. Furthermore, specific abundance ratios constrain various 
nucleosynthetic processes in the young Universe. For example, the abundance of 
the the odd-Z elements P, Ga, and Mn compared to Si, Ge, and Fe, indicates 
an enhanced " odd-even effect" and impact theories of explosive nucleosynthesis. 
Similarly, measurements of the B/O ratio help develop theories of light element 
nucleosynthesis while constraints on Sn, Kr, and other heavy elements will test 
scenarios of the r and s-process. 

4 IONIZED GAS 

The best evidence for ionized gas in damped Lya systems comes from the detec- 
tion of CIV AA 1548.1, 1550.7 in every object for which accurate spectra have 
been obtained. Examples of the C IV velocity profiles are shown in Wolfe & 
Prochaska (2000a) who compare them with the corresponding low-ion profiles. 
The absence of a one-to-one alignment between the velocity components indicates 
that the gas producing C IV absorption is not the same gas that produces low-ion 
absorption. The difference is not surprising. While the X-ray background at z 
~ 3 is about 30 times brighter than at z = (Haardt & Madau 2003), it is still 
not sufficiently bright to produce C 3+ ions by photoionization in neutral gas with 
Lyman limit optical depths, tll >> 10 3 (Wolfe et al. 2004). Note, collisional ion- 
ization is ruled out, as many of the C IV profiles exhibit components with a v < 5 
km s _1 indicating T < 3xl0 4 K, which is lower than the T > 4xl0 4 K threshold 
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required to produce significant fractions of C 3+ . In that respect damped Lya 
systems differ from the ISM in which (a) C 3+ is collisionally ionized in gas with 
T~10 5 K (Sembach & Savage 1992), and (b) the velocity components of the C 
IV and low- ions are aligned. This is interpreted as evidence for corotation of 
hot halo gas with the neutral disk (Savage, Edgar, Sz Diplas 1990). By contrast, 
the C 3+ ions in damped Lya systems are more likely produced in gas that is 
photoionized, has a temperature T~10 4 K, and is kinematically distinct from the 
neutral gas. 

Damped Lya systems also exhibit Si IV AA1393.7, 1402.7 and Al III AA 1854.7, 
1862.7 absorption. Because the Si IV and C IV velocity components are closely 
aligned, the Si +3 ions must be produced by photoionization in the same gas 
containing the C 3+ ions (Wolfe & Prochaska 2000a). In the ISM the Al III lines 
arise in a warm ionized medium (WIM), which is an extensive region of warm 
(T~10 4 K) photoionized hydrogen that pervades the disk of the Galaxy (Reynolds 
2004). Therefore, it is surprising that in damped Lya systems the Al III velocity 
components are aligned with the low ions and misaligned with the C IV and 
Si IV components which arise in gas resembling a WIM. This implies that the 
Al 2+ and low ions both arise either in neutral gas or they arise in photoionized 
gas, which is kinematically distinct from the C 3+ bearing gas. Photoionization 
equilibrium calculations indicate that in most damped Lya systems both Al 2+ 
and the low ions arise in gas that is mainly neutral (Vladilo et al. 2000; Prochaska 
et al. 2002b; Wolfe et al. 2004). The soft X-ray background at z ~ 3 is sufficiently 
bright to produce Al 2+ by photoionization and to leave hydrogen mainly neutral 
and carbon singly ionized (see Wolfe et al. 2004). On the other hand in about 10 
% of damped Lya systems the regions giving rise to low ion and Al III absorption 
contain hydrogen that is significantly ionized. These cases can be recognized by 
the large Fe 2+ /Fe + ratios and low Ar°/Si + ratios (Prochaska et al. 2002b). 

These studies raise some interesting questions. First, photoionization by X- 
ray background radiation produces ions X J for which the ionization potential of 
the preceding ionization state IP(X J_1 ) > 1 Ryd. In many damped Lya systems 
photoionization by locally generated radiation with hv < 1 Ryd produces low ions 
for which IP(X J ' -1 ) < 1 Ryd (see § 8). The question is why doesn't local radiation 
with hv > 1 Ryd photoionize sufficient neutral gas to generate a co-rotating 
WIM like that in the Galaxy? The answer may be related to the higher neutral 
gas content of damped Lya systems, which could reduce the escape fraction 
of ionizing radiation from H II regions in damped Lya systems. Second, do 
damped Lya systems contain hot (T > 10 5 K) collisionally ionized gas? Hot gas 
is a byproduct of feedback processes such as galactic winds or shock heating by 
supernova remnants (Ferrara & Salvaterra 2004). Because of the evidence for 
type II supernovae (see § 3.3.2) and the possible existence of winds (see § 2.5.2), 
hot gas could be present in damped Lya systems. The most efficient technique 
for finding hot gas is through the detection of the O VI AA 1031.9, 1037.6 doublet 
(Sembach et al. 2003). But no surveys for OVI in damped Lya systems have 
been carried out. 

5 MOLECULAR GAS 

Molecular gas is ubiquitous throughout the Galaxy ISM. In a survey for Lyman 
and Werner band absorption, the most efficient tracers of H2 molecules in diffuse 
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interstellar gas, Savage et al. (1977) detected H2 in 90 out of 103 sightlines toward 
background Galactic stars. More specifically, Savage et al. (1977) showed that 
the molecular fraction, /(H 2 ) =2iV(H 2 )/[2iV(H 2 )+iV(H I)], undergoes a steep 
transition from /(H 2 ) < 10~ 4 at 7V(H I) < 4xl0 20 cm" 2 to /(H 2 ) > 10" 2 at 
7V(H I) > 4xl0 20 cm* 2 . Because of the similarity with the 7V(H I)>2xl0 20 cm -2 
threshold, one might expect to find /(H 2 ) > 10~ 2 in a significant fraction of 
damped Lya systems. Furthermore, if damped Lya systems are the neutral 
gas reservoirs for star formation at high redshifts, and since stars form out of 
molecular clouds, molecules should be present. 

However, the H 2 content of damped Lya systems is much lower than in the 
Galaxy. In their compilation of accurate searches for H 2 , Ledoux, Petitjean, & 
Srianand (2003) report the detection of H 2 in only five out of 23 cases qualifying 
as confirmed damped Lya systems, which brings to mind the H 2 content of the 
LMC and SMC. In particular Tumlinson et al. (2002) found H 2 in only 50% of 
the sightlines through the LMC. Moreover, the LMC resembles the damped Lya 
systems in that the typical upper limits are /(H 2 ) < 10 -5 . In addition, the mean 
value of /(H 2 ) for the positive detections is about 10~ 2 for the damped Lya 
systems, the LMC and the SMC, which is about a factor of 10 lower than in the 
Galaxy. 

Why is the H 2 content in damped Lya systems so low? The answer is par- 
tially related to low dust content. In the ISM, H 2 forms on the surfaces of dust 
grains and is destroyed by photodissociation due to FUV radiation. In that case 
/(H 2 )=2i?nn/-f (Jura 1974), where R is the formation rate constant and I is the 
photodissociation rate. Because R oc nn 2 H , /(H 2 ) is predicted to decrease with 
decreasing dust-to-gas ratio, and since / oc J u , where J u is the mean intensity of 
FUV radiation, /(H 2 ) should decrease with increasing radiation intensity. The 
low molecular fractions of the LMC and SMC are plausibly attributed to low 
dust content and high radiation intensities. Similarly, the low dust content of 
damped Lya systems helps to explain the low values of /(H 2 ). Indeed Ledoux, 
Bergeron, & Petitjean (2002) find a statistically significant positive correlation 
between /(H 2 ) and k. 

What is the value of J v in damped Lya systems? Levshakov et al. (2002) 
inferred J u for DLA0347— 38 at z = 3.025 by showing that FUV pumping was re- 
sponsible for populating several excited rotational levels in the ground electronic 
state. By combining the excitation equations with the formation equations it is 
possible to deduce /, independent of the functional form of R, and from that, J v . 
Levshakov et al. (2003) found that J v was comparable to the ambient interstellar 
radiation intensity in the Galaxy, i.e., J^IO -19 ergs cm -2 s _1 Hz -1 sr _1 for this 
system. Previous authors reached similar conclusions for other damped Lya sys- 
tems (see Black, Chaffee, &: Foltz 1987; Petitjean, Srianand, &: Ledoux 2000; Ge 
& Bechtold 1997). This is an important result, because intensities of this mag- 
nitude are significantly higher than the background radiation intensity predicted 
at hv ~ 10 eV and z ~ 3 (Haardt & Madau 2003). Therefore, a local source 
of radiation is required to maintain the right balance between H 2 formation and 
destruction. 

While the molecular content of the diffuse gas detected in damped Lya systems 
is low, dense molecular clouds with high dust content could be present. Such 
objects would be missed owing to obscuration of the background QSOs or due 
to a low covering factor. While future surveys for radio-selected damped Lya 
systems may eventually rule out such scenarios, they are consistent with the 
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current data. In fact dense molecular clouds may be required as the sites of the 
star formation which has been inferred for damped Lya systems. 

6 KINEMATICS 

As described in the previous sections, the introduction of echelle spectrometers 
on 8 to 10 m-class telescopes has led to precise column density measurements 
of weak transitions in the damped Lya systems. It has also led to another 
significant advance in damped Lya research, namely the resolution of kinematic 
characteristics from the velocity profiles of the metal lines. In § 3.1 we showed that 
damped Lya profiles exhibit multiple components that are qualitatively similar 
to the component structure of similar transitions observed in the Galaxy ISM 
(e.g. Wolfe et al. 1994). At high signal-to-noise ratios the line-profiles generally 
decompose into 5 to 30 velocity components, i.e., 'clouds', with column densities 
spanning roughly an order of magnitude. The velocity components comprising 
the low-ion profiles (e.g. Si II A 1808.0) are typically broadened by turbulent 
motions and have velocity dispersions of a v ps 4—7 km s _1 . These characteristics 
are also observed for the high-ion profiles (e.g. C IV A 1548.1; Wolfe & Prochaska 
2000a), although the velocity dispersions are generally larger. 

Prochaska & Wolfe (1997) presented the first modest sample of measurements 
on the low-ion kinematic characteristics of the damped Lya systems. Their results 
and subsequent surveys have shown that the damped Lya systems exhibit velocity 
widths Av ranging from 15 km s^ 1 to several hundred km s _1 with a median of 
ps 90 km s- 1 (Figure 10). Prochaska & Wolfe (1997, 1998) demonstrated that the 
observed Av distribution matched that predicted for rotating disks with typical 
velocity speed v c ~ 200 km s _1 under the important assumption that the gas 
disk is thick (h > O.lRd), where h is the disk scale height and Rd is the radial 
exponential scale length. The authors also stressed that the damped Lya line- 
profiles tend to show the 'edge-leading asymmetry' expected for rotating disks. 
The observations, therefore, suggested a population of disk galaxies with rotation 
speeds corresponding to present-day galaxies. 

The difficulty with this scenario, however, is that hierarchical cosmology im- 
plies that galaxies in the young Universe have smaller masses and lower rotation 
speeds (on average) than the current population. Indeed, Prochaska & Wolfe 
(1997) emphasized that the circular velocity distribution of damped Lya systems 
predicted within the CDM cosmogony (Kauffmann 1996) is incompatible with the 
observations. This includes models (e.g. Mo, Mao, & White 1998) which allow 
for a cross-section weighted distribution of spin parameters (Wolfe & Prochaska 
2000a). As emphasized by Jedamzik & Prochaska (1998) and Prochaska & Wolfe 
(2001), there is a tension within CDM theory between including enough low mass 
galactic halos to match the observed incidence of damped Lya systems without 
severely underpredicting the Av distribution. Granted the many successes of 
CDM theory, perhaps it will remain a true coincidence that the damped Lya 
kinematic characteristics are best described by a population of galaxies similar 
to (albeit thicker than) the present-day disk population. 

If we are to adopt the ACDM power spectrum at z ~ 3, then the velocity 
fields of a significant fraction of damped Lya systems must have contributions 
from non-rotational dynamics in order to explain the kinematic data. Numerical 
simulations within the CDM context describe the damped Lya systems as mul- 
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tiple 'proto-galactic clumps' bound to a virialized dark matter halo (Haehnelt, 
Steinmetz, & Rauch 1998). The kinematics within this scenario are due to the 
combination of infall, random motions, and rotational dynamics. In order for 
this model to be consistent with the observations, however, one requires that 
the damped Lya cross-section A is proportional to v° with a > 2 (Haehnelt, 
Steinmetz, & Rauch 1998; Mailer et al. 2001). As discussed in § 2.5.2, early 
numerical work indicated A oc v^' 1 which implied a crisis between theory and 
observation (Prochaska & Wolfe 2001). More recent results, however, support 
A oc v^- 5 (Nagamine, Springel, & Hernquist 2004a) and the protogalactic clump 
scenario remains a viable option 5 . We stress, however, that no cosmological sim- 
ulation to date has self-consistently matched the low-ion damped Lya kinematic 
observations. 

Lu et al. (1996) first remarked that the high-ion profiles (e.g. CIV A 1548.1) 
of the damped Lya systems have significantly different kinematic characteristics 
from the low-ion gas. Wolfe & Prochaska (2000a) examined a large sample of 
high-ion damped Lya profiles and stressed that while the component structure is 
generally disjoint from the low-ion profiles, the high-ions are roughly centered on 
the low-ion transitions. Furthermore, there is a connection between the velocity 
fields of the two ionization states in that the high-ion gas nearly always shows 
comparable or larger velocity width than the low-ion gas (Figure 11). These 
trends place important constraints on the nature of the damped Lya systems. 
Wolfe & Prochaska (2000b) considered a simple model where the low-ion gas is 
confined to a disk enshrouded in a halo of high-ion gas with kinematics described 
by the infall model of Mo & Miralda-Escude (1996). Wolfe & Prochaska (2000b) 
demonstrated that this model could not simultaneously match the low and high- 
ion kinematics, in particular the co-alignment of the profile centers. In contrast, 
Mailer et al. (2003) demonstrated that a damped Lya model based on multiple 
satellites bound to a single dark matter halo can satisfy the low-ion and high-ion 
kinematics provided each satellite has a halo of C IV producing gas. Their model 
only marginally reproduced the correlation between low-ion and high-ion gas and 
the authors suggested that a kinematic correlation exists beyond their simple 
model. Perhaps numerical simulations of the hot and cold phases will show that 
the hot gas co-rotates with the cold gas. 

Because damped Lya systems are selected solely on the basis of H I column 
density, it is possible that other mechanisms contribute to their velocity fields. 
This could include winds induced by mergers or star formation (Nulsen, Barcons, 
& Fabian 1998; Schaye 2001) or even Hubble flow motions in the ambient IGM. At 
present such scenarios have not been quantitatively developed or tested against 
observation. 

One can gain further insight into the nature of damped Lya systems by synthe- 
sizing the results on kinematic characteristics with their other properties. Wolfe 
& Prochaska (1998) first noted relationships between the velocity width and the 
H I column density and metallicity of the damped Lya systems. They found that 
the damped Lya systems with larger velocity width tend to have higher metallic- 
ity. This trend matches one's physical intuition: If the velocity width is correlated 
with the galactic mass, a correlation between Av and [M/H] is natural provided 
more massive galaxies have higher metallicity. In contrast to the metallicity- A v 

5 We note that in the clump model the 'edge-leading asymmetry' can be reproduced if 3 or 
fewer clumps are intercepted by the majority of damped Lya sightlines. 
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distribution, the damped Lya systems show smaller N(K I) at large Av values. 
The observations run contrary to expectation for disk galaxies where the gradient 
of rotation projected along the line of sight peaks toward the center. The proto- 
galactic clump scenario also would predict a mild, positive correlation (e.g. Mailer 
et al. 2001) or no correlation depending on whether H I surface density correlates 
with halo mass. Therefore, the damped Lya observations are unexpected and dif- 
ficult to interpet in terms of single or multiple rotating disks. Additional studies, 
both observational and theoretical, on correlations between the kinematics and 
other damped Lya properties would place further constraints on the processes of 
galaxy formation in the young Universe. 

It is important to consider the implications of the damped Lya kinematic 
characteristics with respect to the observed relative abundance patterns (e.g. 
Prochaska 2003 ). As noted in § 3.3.2, several authors have recently argued that 
dust-corrected a/Fe ratios of damped Lya systems are roughly solar and indi- 
cate star formation histories representative of dwarf or irregular galaxies instead 
of massive systems (e.g. Calura, Matteucci, k, Vladilo 2003; Tolstoy et al. 2003). 
We note, however, that the velocity widths of the majority of damped Lya sys- 
tems exceed 60 km s _1 and cannot be attributed to the gravitational velocity 
field of a single dwarf galaxy. Furthermore, there is no apparent correlation be- 
tween the gas-phase Si/Fe ratio and velocity width. Therefore, a contradiction 
exists between the observed velocity widths and the interpretation of the damped 
Lya abundance patterns in terms of absorption by a single dwarf galaxy along 
the line of sight. Regarding the CDM clump scenarios, in which multiple dwarf 
galaxies are found along the line of sight, one notes that the majority of damped 
Lya systems are predicted to be embedded in dark matter halos which exceed 
the masses of present-day dwarf galaxies (e.g. Mailer et al. 2001; Nagamine, 
Springel, & Hernquist 2004a) and instead correspond to the progenitors of galax- 
ies like the Milky Way (Steinmetz 2003). It would be particularly valuable to 
perform a simulation which traced the star formation history and resolved the 
galaxy kinematics within a cosmological context to examine abundance pattern 
correlations with gas kinematics. 

7 GALAXY IDENTIFICATIONS 

In this section we describe the results of searches for galaxies physically associated 
with damped Lya systems. 

7.1 Galaxies with z > 1.6 

Because of the presence of bright, B ~ 18.5, background QSOs, surveys for 
galaxies associated with high-z damped Lya systems are more challenging than 
surveys for randomly selected galaxies. To illustrate this point consider a sightline 
passing through an galaxy 6 at a reasonable impact parameter of 10 kpc. At 
z = 3 the galaxy will have an AB magnitude of 24.7 and impact parameter 
of 1.3 arcsec. Detection of the galaxy against the QSO PSF with ground based 
telescopes would be exceedingly difficult even under the best seeing conditions and 
has even proven difficult with broadband space imaging. Using the NICMOS IR 
camera on the HST, Colbert & Malkan (2002) surveyed 22 damped Lya systems 

6 As defined in the z = 3 Lyman Break Galaxy luminosity function of Steidel et al. (1999). 
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and detected only one candidate counterpart down to Hab = 23.5, implying that 
most damped Lya systems are not drawn from the luminous end of the Lyman 
Break Galaxy luminosity function. Warren et al. (2001) probed even deeper, to 
Hab = 25, and found 41 candidate counterparts near 18 high-redshift damped 
Lya systems. Broadband imaging to the necessary depth is thus limited by source 
confusion within reasonable impact parameters. 

The most widely used techniques are searches for Lya emission lines at the 
absorption rcdshift. The advantage of this method is that the wavelength of Lya 
emission is located at the bottom of the damped Lya absorption trough, which 
blocks the bright light of the background QSO. As a result, background night 
sky emission is the only source of external noise. Using slit spectroscopy, Foltz, 
Chaffee, &: Weymann (1986) failed to detect Lya emission with a 3-ct upper 
limit of F ~ 10~ 16 ergs cm -2 s _1 for an unresolved object. While Hunstead, 
Fletcher, & Pettini (1990) claimed detection of Lya emission from a compact 
source coinciding with DLA0836+11 at z = 2.466, this feature was not confirmed 
in spectra acquired by Wolfe et al. (1992) and Lowenthal et al. (1995). Nor was 
Lya emission detected in imaging surveys using narrow-band interference filters 
or Fabry-Perot interferometers. In this case the QSO light is blocked because the 
bandwidth of the filter is centered on the damped Lya line but has a narrower 
FWHM. This technique is ideal for impact parameters large compared to the 
seeing radius, since an emitter located outside a slit could still be detected in the 
narrow band image. Smith et al. (1989), Deharveng, Bowyer, & Buat (1990) and 
Wolfe et al. (1992) carried out narrow-band surveys for Lya emission. Lowenthal 
et al. (1995). carried out Fabry-Perot surveys. No detections to limiting fluxes of 
F~5xW~ ergs cm~ 2 s _1 for unresolved objects were reported. The extended 
Lya emitter associated with DLA0836+11 at z = 2.466 claimed by Wolfe et 
al. (1992) is more likely to be a galaxy associated with a lower redshift Mg II 
absorption system (Lowenthal et al. 1995). 

The failure to detect Lya emission could result from the destruction of reso- 
nantly trapped photons by even a small amount of dust (Chariot & Fall 1991). 
For this reason several groups attempted to detect damped Lya systems in Ha 
emission since Ha photons are not resonantly trapped. Bunker et al. (1999) 
used IR detectors on a 4 m class ground-based telescope to search for Ha in five 
damped Lya systems, but none were detected (see also Mannucci et al. 1998). 
Of course, the failure to detect Ha emission might also be due to small impact 
parameters, since an emitting region located within the PSF of the QSO would 
not be detected from the ground. However, Kulkarni et al. (2000, 2001) used 
NICMOS to search for Ha emission from two damped Lya systems and none was 
found (see Table 1). 

With the recent detections of Lya emission from at least 2 out of 18 damped 
Lya systems surveyed using 8-to 10-m class telescopes (M0ller et al. 2002; M0ller, 
Fynbo, Sz Fall 2004), it is evident that the previous null detections of Lya were 
largely due to the lower sensitivity of 4-m class telescopes. The results, summa- 
rized in Table 1, indicate that two of the three known Lya emitters, DLA0458— 02 
and DLA0953+47A, would not have been detected in the earlier surveys. While 
the sample is too small to draw general conclusions, the results are interesting 
for the following reasons. First, damped Lya systems resemble randomly se- 
lected Lya emitters by the similarity in Lya luminosity and in the compact size 
of the emission regions. Second, the small impact parameters for DLA0458— 02 
and DLA0953+47A suggest that the H I absorbing layers are smaller than ~ 5 
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kpc. However, since the H I content of DLA0458— 02 is known to extend over 
linear scales exceeding 17 kpc (Briggs et al. 1989), the sightline to the QSO must 
pass close to a compact star forming region, which is embedded in a much larger 
layer of H I. Third, no continuum emission has been detected from the same two 
damped Lya systems. The B > 27 limit on DLA0953+47A places this damped 
Lya system near the faint end of the known luminosity function of Lyman Break 
galaxies. 

By contrast, DLA2206— 19A is a luminous Lyman Break galaxy. This emitter 
was first located in IR images obtained with NICMOS (Warren et al. 2001). The 
STIS image (M0ller et al. 2002) shows rest-frame FUV stellar emission extending 
1 arcsec between a bright knot and the QSO sightline. Lya emission at the 
redshift of the damped Lya system was detected from the knot with spectra 
obtained with the VLT (M0ller et al. 2002). The magnitude integrated over the 
object is V= 23 (M0ller 2004, priv. comm.), which places this damped Lya 
system at the bright end of the Lyman Break luminosity function. 

As a result, there is little overlap between the luminosity functions of damped 
Lya systems and the R < 25.5 spectroscopic sample of Lyman Break galaxies (see 
M0ller et al. 2002; Schaye 2001) Efforts to detect the clustering of damped Lya 
systems with neighboring Lyman Break galaxies have so far only yielded upper 
limits (Gawiser et al. 2001; Adelberger et al. 2003), providing further evidence 
that the damped Lya systems are nearly disjoint from the R < 25.5 spectroscopic 
sample of Lyman Break Galaxies, with DLA2206— 19A a clear exception. Based 
on the UV continua implied by the Lya luminosities of the other 2 damped 
Lya systems detected in emission, there is growing evidence that at least some 
damped Lya systems overlap with the dimmer "photometric" sample of Lyman 
Break Galaxies at R < 27 seen in the Hubble Deep Fields. 

7.2 Galaxies with z < 1.6 

While the nature of the galaxies associated with damped Lya systems at z > 1.6 
is still unclear, at lower redshifts there should be a close resemblance to objects 
drawn from the population of normal galaxies if damped Lya systems trace the 
star formation history of normal galaxies. However, the low metallicities inferred 
for most low-z damped Lya systems (Pettini et al. 1999; Kulkarni et al. 2005) 
has cast doubt on this idea and has led to the suggestion that low-z damped 
Lya systems are metal-poor objects such as dwarf galaxies (Calura, Matteucci, 
& Vladilo 2003) or low surface-brightness galaxies (Jimenez, Bowen, & Matteucci 
1999). As a result, identification of galaxies associated with damped Lya systems 
at z < 1 is of vital importance. 

There are 23 damped Lya systems known at z < 1.6 (as of October 2004). 
Thirteen of these have been identified using spectroscopic or photometric red- 
shift techniques (Chen & Lanzetta 2003; Chen, Kennicutt, & Rauch 2005) and 
possible host galaxies for an additional eight systems have been found in imag- 
ing surveys (Le Brun et al. 1997; Rao et al. 2003). Chen & Lanzetta (2003) 
analyzed an unbiased subset of nine galaxies with redshifts and concluded that 
the resultant luminosity distribution was dominated by luminous galaxies with 
L/L^y 0.1; i.e., luminous galaxies dominate the neutral gas cross-section at z < 
1. More specifically, they used a maximum- likelihood technique to determine the 
dependence of H I cross section on luminosity and showed that a cross-section 
weighted Schechter function with typical parameters for normal galaxies provided 
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a good fit to the data. 

On the other hand, Rao et al. (2003) analyzed a heterogeneous sample of 
fourteen damped Lya system host galaxy candidates, including objects without 
confirmed redshifts and concluded that the neutral gas cross-section at z < 1 
was dominated by dwarf galaxies. Since Rao et al. (2003) did not quantify their 
result, this conclusion is difficult to evaluate. However, comparison with the 
sample of Chen & Lanzetta (2003) shows a larger fraction of galaxies with lower 
values of L/L* in the Rao et al. (2003) sample. In some cases these are galaxies 
without confirmed redshifts at small angular separations from the damped Lya 
sightline. In other instances, where two or more galaxies are found to be associ- 
ated with the damped Lya system, Rao et al. (2003) chose the lowest luminosity 
galaxy for the analysis because it had the smallest impact parameter. However, 
Chen & Lanzetta (2003) argue that in these cases, an entire group of galaxies is 
responsible for the absorption profile, so the object with lowest impact parame- 
ter may not be the appropriate choice. Given these uncertainties and the small 
sizes of the samples, we conclude that the current data are consistent with the 
galaxies associated with low redshift damped Lya systems being drawn from the 
population of normal galaxies. 

8 STAR FORMATION IN DAMPED Lya SYSTEMS 

Several lines of evidence imply that damped Lya systems experience ongoing 
star formation. At z < 1.6 the evidence is unambiguous, since more than half the 
sample is associated with galaxies of stars. At z > 1.6 the evidence is indirect 
because only one damped Lya system is identified with a resolved object of 
galactic dimensions that emits starlight. While starlight likely ionizes the gas 
that gives rise to Lya emission in two other high-z damped Lya systems, the 
starlight itself has not been directly detected. The expectation is that most of 
the metals seen in damped Lya systems were generated through star formation 
in their host systems, but the instantaneous star formation rates characterizing 
this important population of objects remained unknown until recently. 

8.1 Direct Emission Measurements of Star Formation Rates in Damped 
Lya Systems with z > 1.6 

Table 1 summarizes results including SFRs for the three positive detections and 
for objects with upper limits on SFR at z > 1.6. Observed lower limits on L(Ha) 
were coverted to upper limits on SFR using the Kennicutt (1998) calibration, 
SFR(M yr- 1 )=7.9xlO~ 42 L(Ha)(ergs s" 1 ). Detections of Lya emission 7 were 
combined with the expression for L(Ha) and case B radiative recombination to 
find lower limits of SFR(M Q yr _1 )=l.lx 10 _42 L(Lya)(ergs s _1 ). The total magni- 
tude integrated over DLA2206— 19A of V= 23 and the very sensitive B > 27 limit 
on DLA0953+47A (A. Bunker 2004, priv. comm.) were combined with the rest- 
frame UV Kennicutt (1998) calibration, SFR(M yr _1 )=1.4x 10^ 28 L J/ (ergs s" 1 Hz" 1 ), 
to set a lower limit on the SFR for DLA2206— 19A and a rough upper limit for 
DLA0953+47A assuming no dust extinction for the latter object. An estimated 
dust extinction correction for DLA2206— 19A (Wolfe et al. 2004), suggests an 

7 Due to its extreme sensitivity to dust extinction, Lya emission provides a lower limit to the 
star formation rate as long as the AGN contribution to the emission line is negligible. 
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upper limit on SFR of 50 M yr 1 , which is comparable to the SFRs of the more 
luminous Lyman Break galaxies (Shapley et al. 2003). 

8.2 Star Formation Rates from the C II* Technique 

A method recently developed by Wolfe, Prochaska, & Gawiser (2003b) makes it 
possible to infer the star formation rate per unit area for individual damped Lya 
systems. The basic method is to use measurements of the C II* 1335.7 column 
density to measure the [C II] 158 um cooling rate in the neutral gas producing the 
damped Lya absorption. This is possible because the C II* A 1335.7 transition 
arises from the excited 2 Pz/2 state in C + , and spontaneous photon decay of the 
2 -P 3 / 2 state to the 2 P\/2 state results in [C II] 158 fim emission, which is the 
principal coolant of neutral gas in the Galactic ISM (Wright et al. 1991). Under 
the presumed condition of thermal balance, the cooling rate equals the heating 
rate and it is possible to calculate the star formation rate per unit area that 
generates the implied heating rate. 

By analogy with the Galactic ISM, Wolfe, Prochaska, & Gawiser (2003b) adopt 
the grain photoelectric effect as the principal heating mechanism for damped Lya 
systems. In that case FUV radiation {hu^, 6 to 13.6 eV) ejects photoelectrons 
from grain surfaces, which heat ambient electrons through Coulomb interactions 
(Bakes &: Tielens 1994; Weingartner & Draine 2003). The heating rate from 
the grain photoelectric effect is proportional to the FUV radiation intensity J u , 
which consists of a contribution from the FUV background radiation plus a local 
contribution from hot stars located inside the galaxy which is proportional to the 
instantaneous star formation rate per unit area, ip*. 

To infer ip* and other properties from observations, Wolfe, Prochaska, & Ga- 
wiser (2003b) deduce the [C II] 158 um spontaneous emission rate per H atom 
from the observational quantity: 

. _ N(C U*)hu ul A ul _! 



7V(H I) 



ergss^BT 1 , (8) 



where A u \ is the Einstein coefficient and hv u \ is the energy of the 158 ^m transi- 
tion. These are known quantities, and N(C II*) and iV(B I) are measured from 
the absorption line spectra. The total heating rate includes inputs due to the 
grain photoelectric effect, X-ray photoionization, cosmic ray ionization, C I pho- 
toionization, and collisional heating. The total cooling rate includes cooling due 
to Lya, grain radiative recombination, and emission by fine-structure states of 
O , Si + , Fe + , etc. in addition to the [C II] 158 //m fine-structure emission. Ra- 
diative excitation of the C + fine-structure states by CMB radiation is included 
because it can be significant at high redshift. Wolfe et al. (2004) showed that 
the heating rates predicted for the Baardt & Madau (2003) backgrounds were 
significantly lower than the 158 um. cooling rates implied for damped Lya sys- 
tems with detected C II* A 1335.7 absorption, requiring active star formation 
to explain the observed cooling rates. On the other hand, Wolfe et al. (2004) 
also showed that heating by background radiation alone cannot be ruled out for 
systems in which C II* A 1335.7 absorption was not detected. A summary of the 
l c data for damped Lya systems (and the Galaxy) is given in Figure 12. Wolfe, 
Prochaska, & Gawiser (2003b) argued that the heating rate in the Galaxy is much 
larger than in damped Lya systems because the dust content in the Galaxy is 
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at least a factor of 30 higher than in damped Lya systems, whereas V* i n the 
Galaxy is only a factor of 2 to 3 lower than in damped Lya systems (see below). 

In the case of systems with detected C II* A 1335.7 absorption, one compares 
158 /xm emission rates computed for a range of star formation rates per unit 
area with the empirical quantity, £ c . This results in in two solutions for ip* 
corresponding to thermally stable states of a two-phase medium in which warm 
neutral medium (WNM) gas is in pressure equilibrium with cold neutral medium 
(CNM) gas: the WNM solution corresponds to gas with T~8000K and n~0.02 
cm~ 3 , and the CNM solution corresponds to gas with T~100K and n~10cm~ 3 . 
In every case the WNM solution for ip* is a factor of ten or more higher than 
the CNM solution: since C II emission in the WNM is a small fraction of the 
total cooling rate, the total heating rate implied for an observed value of i c must 
be higher in the WNM than in the CNM. Owing to the low dust optical depths 
of damped Lya systems the resulting tp* are system-wide average star formation 
rates rather than local values as is the case in the dusty Galactic ISM. 

This method was used by Wolfe et al. (2004) to analyze a sample of 45 damped 
Lya systems, 23 with measured C II* column densities and 22 for which C II* 
was not detected. In the redshift interval z =[1.6,4.5] the average SFR per unit 
area for positive detections is <^*> = 11.3 xlO -3 M yr _1 kpc -2 for the CNM 
solution and 0.21 M yr _1 kpc -2 for the WNM solution; by comparison in the 
Galaxy tjj* ~ 4xl0 -3 Mq yr _1 kpc -2 . However, the WNM solution is unlikely 
to be correct, because the bolometric background intensity produced exceeds the 
observational limits (see Wolfe, Gawiser, Sz Prochaska 2003a) . This conclu- 
sion is consistent with recent evidence from individual systems; Howk, Wolfe, & 
Prochaska (2005) showed that the low upper limit on the optical-depth ratio of 
Si II* A 1264.7 to C II* A 1335.7 in a high-redshift damped Lya system results 
in an upper limit of 800 K for the temperature in the gas producing C II*, which 
implies that C II* absorption arises in CNM gas in this system. Furthermore, a 
pure WNM solution would significantly overpredict the observed rest-frame-UV 
luminosity from DLA2206— 19A based on its measured £ c , implying that C II* 
absorption in this damped Lya system arises in a CNM. It is significant that the 
predicted value for CNM, J v cl1 * = 1.7l 2 ;gXlO" 18 ergs cm~ 2 s _1 Hz -1 sr _1 is the 
largest mean intensity inferred from the entire C II* sample. This may help to 
explain why DLA2206— 19A is one of the rare damped Lya systems detected in 
emission. 

On the other hand the gas detected in absorption is more likely to be a WNM 
for damped Lya systems with upper limits on C II* absorption. In this case the 
gas could be a pure single-phase WNM heated by background radiation alone. 
These damped Lya systems would then be objects without significant star for- 
mation at the absorption epochs. Or the gas could be the WNM branch of a 
two-phase medium in which the SFR per unit area is similar to that found for 
the CNM solutions (Wolfe et al. 2004). While the absence of 21 cm absorption 
at z > 3 (Kanekar & Chengalur 2003) and the large C II/C I ratios (Liszt 2002) 
are consistent with the WNM hypothesis, both phenomena are also naturally ex- 
plained within the context of the two-phase hypothesis. Clearly, it is important 
to determine which of these explanations is the correct one. 
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8. 3 Implications 

The C II* technique can also be used to obtain quantities with cosmological 
significance. Specifically, the SFR per unit comoving volume is given by p\ = 
(tp*(z))(A :¥ /Ap)(dN/dX)( y Ho/c), where A* is the intrinsic cross-sectional area 
occupied by stars emitting FUV radiation and A p is the projection of the intrinsic 
H I area, A, on the sky. 8 In the uniform disk model, neutral gas and stars are 
uniformly distributed across A, yielding the values for p\(z) given in Table 2, 
which also summarizes the systematic uncertainties discussed at length in Wolfe, 
Gawiser, & Prochaska (2003a). 

The total mass of stars produced in damped Lya systems as a function of 
redshift can be computed by integrating over the cosmic star formation history 
of damped Lya systems. The mass density of stars predicted by today, including 
the contribution from "normal" galaxies at z < 1.6, is consistent with the total 
current density of stars in spiral and elliptical galaxies, but this does not reveal 
the precise population of stars produced by high-redshift damped Lya systems. 
The same star formation history will also deplete the neutral gas reservoir at high 
redshift in about ~ 2 Gyr; i.e., £l g {z) would vanish by z ~ 2 if star formation 
starts at z ~ 5. Instead £l g {z) drops by a factor of two during this time, which 
argues for replenishment of neutral gas at an accretion rate p a (z)=0.5p*(z). 

Wolfe, Gawiser, & Prochaska (2003a) also calculated the mass density of met- 
als produced in damped Lya systems by using the conversion of star formation 
rate to metal formation rate suggested by Madau et al. (1996) and Pettini (1999), 
Pmetais = (l/42)p*. Extrapolating to the present day, including the contribution 
from star formation in normal galaxies at z < 1.6, yields a mass of metals a few 
times larger than that found in spiral bulges today, which seems feasible. Inte- 
grating under the DLA cosmic star formation history at z > 2.5 predicts a factor 
of 30 higher cosmic mean metallicity due to metal enrichment of neutral gas than 
is observed in damped Lya systems at z = 2.5 (Prochaska et al. 2003a). This 
sort of "missing metals" problem was first identified by Pettini (1999) for Lyman 
Break Galaxy star formation rates compared to damped Lya system metallicities 
but that problem could be solved by assuming that damped Lya systems are not 
the descendants of the star-forming Lyman Break Galaxies. The damped Lya 
system "missing metals" problem appears to be fundamental and may illumi- 
nate a fundamental flaw in our understanding of the relationship between metal 
enrichment and star formation at high redshift. The problem is a significant 
challenge not only to the C II* method of measuring damped Lya star formation 
rates, but to most hierarchical models, which also produce too many metals (see 
Nagamine, Springel, & Hernquist 2004b; Somerville, Primack, & Faber 2001) 9 . 
Given the other successes of the method it seems likely that the problem has a 
physical solution. 

The simplest solution would be to hypothesize that the damped Lya systems 
forming stars at z > 2.5 are no longer damped Lya systems at z = 2.5 but have 
used up most of their neutral gas and now have the expected metallicities; the 
Lyman Break galaxies are an obvious candidate for their descendants. This is 
inconsistent, however, with timescale of ~ 2 Gyr to use up most of the neutral 

8 Note that A p equals A(N, X) (see eq. 1 ) averaged over the column-density interval JV(H 

9 The model of Pei, Fall, & Hauscr (1999) avoids this problem by invoking significant obscu- 
ration corrections and low yields. 
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gas implied by the observed star formation rates, which is too long. Ejecting the 
metals into the IGM does not solve the problem, as the IGM metallicity predicted 
at z = 2.5 would be [M/H]= —1.5 i.e. about a factor of 30 higher than observed 
in the Lyman a forest, unless the ejected gas is so hot that most metals are in 
an unobservable ionization state, with the possible exception of oxygen, which 
might be detected in O VI AA1031.9, 1037.6. An alternative possibility is that 
the metals are sequestered, with the most attractive solution being a model in 
which damped Lya systems represent neutral gas on the outskirts of actively star- 
forming "bulge" regions with the majority of the produced metals remaining in 
these regions rather than polluting the damped Lya gas or general IGM (Wolfe, 
Gawiser, & Prochaska 2003a). This hypothesis necessitates the conclusion that 
metal-rich bulge regions are underrepresented (or completely unrepresented) in 
the HI- weighted cosmic mean metallicity measured from DLAs. That would not 
be surprising given that the bulges may well have used up most of their neutral 
gas and/or have sufficient dust content to dim background QSOs out of optically 
selected samples. 

9 CHEMICAL EVOLUTION MODELS for DAMPED Lya SYS- 
TEMS 

In the "global" approximation, one averages all quantities over large comoving 
volumes and then solves the chemical evolution equations to deduce metal pro- 
duction rates from the comoving SFR, p*(z), which is computed by tracking 
changes in the neutral gas density, p g {z) (Lanzetta et al. 1995; Pei & Fall 1995; 
Malaney Sz Chaboyer 1996; Edmunds & Phillipps 1997). However, the factor-of- 
two decrease between z ~ 4 and z ~ 2 suggests that damped Lya systems are 
replenished by a net inflow of neutral gas (see § 8.3). 

Pei & Fall (1995) considered models with inflow and outflow and searched for 
self-consistent evolution of the neutral gas, metallicity, and dust. These authors 
calculated the effects of obscuration on quantities deduced directly from the data, 
such as N(R I), and found that significant obscuration was necessary to explain 
the observations. Pei & Fall (1995) fitted the available data with analytic func- 
tions for p g (z) that increased with z at z < 2 (Lanzetta, Wolfe, & Turnshek 1995) 
and assumed that the net accretion rate, p a (z), was proportional to p*(z). When 
more accurate measurements of p*(z) (Steidel et al. 1999) and the cosmic back- 
ground radiation intensity (Hauser & Dwek 2001) became available, Pei, Fall, 
& Hauser (1999) also included the production of background radiation by stars 
in damped Lya systems (see Fall, Chariot, & Pei 1996). With these additional 
constraints Pei, Fall, & Hauser (1999) worked directly from the measurements of 
pg(z) and eliminated the assumption that p a {z) was proportional to p\{z). The 
newer models reproduced the more accurate metallicity measurements not avail- 
able earlier, and were consistent with measurements of the background radiation 
intensities and p*(z). However, as we discuss in § 10, these models appear to 
cause more obscuration than the current observational data imply. Furthermore, 
the values of p\{z) at z <2 were inferred from changes in p g (z) that now appear 
to be spurious (Rao, Turnshek, & Nestor 2004; priv. comm.). 

In the "local" approximation, one computes the chemical evolution of isolated 
galaxies outside a cosmological setting. A star-formation history is imposed from 
the outset and one solves for the chemical response of stars and gas. Adopt- 
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ing the slow star formation history for spiral galaxies suggested by Mateucci, 
Molaro, & Vladilo (1997), Lindner, Fritze-v. Alvensleben, & Fricke (1999) re- 
produced the slow evolution of [Zn/H] with z observed in damped Lya systems. 
However, these authors adopted a "closed box" model and neglected spatial gra- 
dients in all physical parameters. Calura, Matteucci, & Vladilo (2003) removed 
these restrictions and also reproduced the slow increase of [Zn/H] with decreas- 
ing redshift using star formation histories predicted for large galactic radii in 
spirals or for episodic bursts in dwarf irregulars. Furthermore, they used the 
same models to reproduce the [Si/Fe] versus [Fe/H] relation after correcting for 
depletion. Dessauges-Zavadsky et al. (2004) used these models to explain the 
chemical evolution of three damped Lya systems for which abundances of a large 
number of elements had been obtained. Interestingly, the predicted SFRs per 
unit area agree with those inferred from the C II* technique. Whereas Wolfe, 
Gawiser, & Prochaska (2003a) found that integrating such SFRs between z = 5 
and 3 resulted in the overproduction of metals, Dessauges-Zavadsky et al. (2004) 
found that the cumulative metals produced did not exceed those observed owing 
to the short time-scales for metal production required to explain relative abun- 
dance ratios such as [Si/Fe]. However, in some cases the short time scales conflict 
with the conservative lower limit of 0.25 Gyr on age set by the measurement of 
[N/a] near the —0.7 plateau (see § 3.3.2). Other potential problems with these 
models stem from the depletion corrections applied to the [Si/Fe] ratio, which 
may be too large (see § 10). 

The most promising approach to chemical evolution is the direct one, which 
uses cosmological hydrodynamical simulations (see Somerville, Primack, & Faber 
2001 and Mathlin et al. 2001 for semi-analytic and analytic variants of this 
method). The simulations unite the "local" and "global" approximations with 
a self-consistent evolution of stars, gas, metals, and dust within a ACDM cos- 
mology. While the microphysics behind star formation and metal production 
cannot be included in these simulations, recipes calibrated to local observations 
are used to track stars and metals along with dark matter particles governed by 
gravity and gas particles governed by gravity and hydrodynamics. As a result, 
processes such as accretion of neutral gas from the IGM are described physically, 
and star formation is treated self-consistently rather than being imposed ad hoc. 
Moreover, merging between dark-matter halos is included for the first time. 

Cen & Ostriker (1999) were the first to describe the chemical evolution of 
damped Lya systems with numerical simulations. Using a low-resolution Eulerian 
scheme, these authors were unable to resolve the dark-matter halos giving rise to 
damped Lya absorption. Nevertheless, they pointed out that metallicity is a more 
sensitive function of overdensity, 8, than of age: metal-poor objects such as the 
Lya-forest clouds formed in low density environments (S ~ 1), while more metal- 
rich objects such as damped Lya systems and Lyman Break Galaxies formed 
in regions of higher overdensity (5 > 10). Using a more accurate version of 
this numerical code, Cen et al. (2003) predicted the cosmic metallicity at z ~ 3 
to be between 0.3 and 0.5 dex higher than the observed value. They solved 
this "missing metals problem" (see § 8.3) by using obscuration corrections that 
may be larger than allowed by the results of Murphy & Liske (2004). They 
also predicted that, independent of metallicity, the ages of typical damped Lya 
systems in the redshift interval z=[2,4] would be 0.8—2 Gyr, which are consistent 
with the presence of the upper [N/a] plateau. Another prediction of interest is 
that the median stellar mass M*~10 9 M Q , which is a factor of 10 lower than that 
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of Lyman Break galaxies, indicating they are different populations. 

The present state-of-the-art in cosmological hydrodynamic simulations of damped 
Lya systems is represented by the recent results of Nagamine, Springel, & Hern- 
quist (2004b) who used the SPH code described in § 2.5.2 (see also Cora et al. 
2003). These authors found SFRs per unit area that agree with the predictions 
of the C II* technique for damped Lya systems (Wolfe, Prochaska, Sz Gawiser 
2003b). They also predicted an overproduction of metals by z~2.5, but in this 
case by a factor of 10 compared to the observed metal abundances. The difference 
with Cen et al. (2003) is likely related to the lower spatial resolution of the latter 
simulation (about 30/i _1 kpc comoving), which causes the high-metallicities of 
compact regions to be diluted by the low metallicities of diffuse regions. One 
of the interesting predictions of the Nagamine, Springel, & Hernquist (2004b) 
simulations is that all regions in which N(R I) > 2xl0 20 cm -2 exhibit star for- 
mation. Confirmation of this prediction would favor the uniform disk model over 
the bulge model of star formation discussed by Wolfe, Gawiser, & Prochaska 
(2003a). As a result, it is important to decide whether this finding is an artifact 
of the star formation algorithm employed by Nagamine, Springel, & Hernquist 
(2004b), especially since there are regions in nearby galaxies in which iV(H I) 
> 2xl0 20 cm -2 , but only low star formation rates (~10~ 5 M yr -1 kpc -2 ) are 
observed Ferguson et al. (1998). 

10 ARE DAMPED Lya SAMPLES BIASED BY DUST? 

Surveys for damped Lya systems have the greatest impact if they represent a fair 
sample of the neutral gas in the Universe, allowing a clear probe of the evolution 
with redshift of the neutral hydrogen content and the metallicity of neutral gas. 
However, it has long been a major concern that the sample of damped Lya 
systems suffers from "dust bias" i.e. the absence from a magnitude-limited QSO 
sample of those QSOs which suffer obscuration from dusty foreground damped 
Lya systems, leading to underrepresentation of dusty damped Lya systems in the 
overall sample. The easiest way to probe the existence and abundance of dust in 
damped Lya systems would be to find the 2175A bump feature superimposed in 
absorption on background QSO spectra. While at least one strong example has 
been found (Junkkarinen et al. 2004), this does not appear to be the rule (Pei, 
Fall, Sz Bechtold 1991). Without such sharp features to look for and given the 
wide range of intrinsic QSO spectral slopes, reddening from dust in damped Lya 
systems must be searched for statistically by checking if the sample of QSOs with 
foreground damped Lya systems is redder on average than a "control sample" of 
QSOs without foreground damped Lya systems. 

10.1 Observational estimates of reddening 

Ostriker & Heisler (1984) pointed out that optically-selected QSO samples are 
biased towards those QSOs with little foreground dust extinction. Fall & Pei 
(1989) showed that dust in damped Lya systems did not appear to cause the fa- 
mous drop in the QSO number abundance at z > 3. Pei, Fall, & Bechtold (1991) 
detected reddening from damped Lya systems at the \a confidence level and 
inferred dust-to-gas ratios between 1/20 and 1/5 that of the Galaxy, enough to 
explain the lack of observed Lyman a emission from damped Lya systems. This 
led to the prediction that 10%- 70% of QSOs are missing from optically-selected 
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samples, leading to an order of magnitude uncertainty in Q g , < Z >, and other 
quantities estimated from damped Lya systems (Fall & Pei 1993). However, the 
dust-to-gas ratios estimated from high-resolution echelle spectroscopy of QSOs 
with foreground damped Lya systems are lower than the dust-to-gas ratios pre- 
dicted by Fall & Pei (1993), reducing the uncertainty in quantities such as Q g to 
factors of 2-3. Pettini et al. (1997a) combined a metallicity of 1/15 solar with a 
dust-to-metals ratio of 1 /2 that in the Milky Way to find a typical damped Lya 
system dust-to-gas ratio of 1/30 Galactic. Using an SMC reddening curve, they 
predicted a dust extinction of only 0.1 magnitudes at 1500A in the spectrum of 
background QSOs due to damped Lya system dust. If a nucleosynthetic floor 
exists in damped Lya systems at [Si/Fe]~ 0.3 then the dust-to-gas ratios are even 
lower than this, closer to 1/200 Galactic in most systems. Indeed, the detection 
of reddening due to damped Lya systems by Pei, Fall, & Bechtold (1991) con- 
flicts with the recent finding by Murphy & Liske (2004) that E(B -V) < 0.01 
magnitudes using 81 damped Lya systems found in a homogeneous set of SDSS 
Data Release 2 QSOs. The resolution of the conflict is not clear at present. 

10.2 Surveys of radio-selected QSOs 

An insidious (but not physically motivated) possibility would be the existence of 
gray dust associated with the damped Lya systems which could cause obscura- 
tion without the telltale effect of reddening. Even this effect could be overcome 
by using a radio-selected sample of QSOs. The reason this has not typically been 
done is two-fold: (1) The ability to select QSOs within a preferred redshift range 
makes optical color selection more efficient. (2) Conducting optical spectroscopic 
follow-up on a radio-selected sample of QSOs is far more time-consuming pre- 
cisely because they do not have a strict optical magnitude limit; half of the total 
exposure time can be required by the dimmest one or two objects. Obviously, 
dropping those from the survey would defeat the entire purpose of radio selection. 
One radio-selected sample has been published: the CORALS survey of Ellison 
et al. (2001) found 19 intervening damped Lya systems towards 66 z em > 2.2 
radio-selected QSOs from the Parkes quarter- Jansky sample (Shaver et al. 1996), 
yielding a marginal increase in dN /dX and Cl g at < z >= 2.37 versus optically- 
selected QSO samples. These results imply that at most half of damped Lya 
systems are missing from optically-selected samples. For Q g , the radio sample 
yields 1.4xl0~ 3 as opposed to the value of 6.7xl0 -4 found for optically selected 
samples at this redshift (Prochaska & Herbert-Fort 2004), but this is only a 1.5 
a difference given the small sample size. 

10.3 Empirical Estimates of Damped Lya System Obscuration 

A third way to estimate the effects of dust obscuration by damped Lya systems 
is to infer this from the observed chemical abundances. Taking the observed HI 
column densities and the dust-to-gas ratios implied by the depletion patterns of 
the damped Lya systems (see Eq. 7), it is possible to estimate the extinction in 
the rest-frame UV of the QSO for an assumed extinction curve. Given the lack of 
the observed 2175A bump feature, it appears more reasonable to assume an SMC 
(Prevot et al. 1984) rather than Galactic (Cardelli, Clayton, & Mathis 1989) dust 
extinction law. Prochaska & Wolfe (2002) used this technique (see their Figure 24; 
see Prochaska 2004 for an update) to correct the observed QSO magnitudes by 
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this inferred extinction and then to compare the implied true magnitude with the 
magnitude limit of the survey used to search for damped Lya systems (which is 
typically shallower than the limit of the survey used to discover the QSOs). These 
quantities were then compared to a bootstrap prediction of how many QSOs are 
expected to be missing from observed samples due to extinction by foreground 
damped Lya systems. The typical range of extinction corrections runs from 
to 0.3 magnitudes, even though half of the QSOs are so much brighter than the 
survey limit that they could have been seen with up to 1 magnitude of extinction. 
This shows that damped Lya systems which cause between 0.3 magnitudes and 

I magnitude of extinction are rare and predicts that at most 10% of QSOs are 
missing from the samples probed for damped Lya systems due to a damped Lya 
dust bias. 

II ARE DAMPED Lya SAMPLES BIASED BY GRAVITA- 
TIONAL LENSING ? 

On the other hand Prochaska, Herbert-Fort, & Wolfe (2005) find that biasing 
due to gravitational lensing could be important. They compared the full SDSS 
sample of damped Lya systems with subsamples comprising the brightest 33 % 
of background quasars and the faintest 33 % of background quasars. While the 
incidence of damped Lya systems, dM /dX was found to be insensitive to quasar 
magnitude, the mass density, Cl g (z), was found to vary significantly. Specifically, 
the bright subsample showed systematically higher values of £l g (z) than the faint 
subsample. To explain the independence of dM jdX on quasar magnitude, the 
difference must lie in the incidence of systems with large iV(H I), which is observed 
to be larger in the bright subsample. 

Prochaska, Herbert-Fort, & Wolfe (2005) argue that gravitational magnifica- 
tion of the background quasars by massive halos or disks associated the fore- 
ground damped Lya systems could account for this systematic effect, which was 
first detected at the 2a level by Murphy & Liske (2004). Whereas obscuration 
by dust would cause £l g (z) to be lower in the bright subsample, magnification by 
lensing due to exponential disks is greatest for damped Lya systems with large 
values of N(K I) (Bartelmann & Loeb 1996; Mailer et al. 1997). 

Consequently, the values of Q g (z) in Figure (5) may be 10-20 % too high, but 
the evolution of Q g (z) with z is still likely to be correct. However, the effect could 
have a more important effect on damped Lya system sample with z <2 (Rao & 
Turnshek 2000). 

12 CONCLUSIONS 

We end this review with the following question: What have we learned from 
damped Lya systems that we did not know before? We attempt to answer this 
question by listing results judged to be robust. These are also summarized in 
Tables 3 and 2, which describe cosmological and local properties, respectively. 
Table 3 lists both the medians and the means in order to show the effects of the 
upper limits placed on various parameters. Specifically, the means only include 
positively detected quantities while the medians include upper limits. 

(1) Most of the neutral gas in the Universe in the redshift interval z= [0,5] 
is in damped Lya systems. The cosmology and mean intensity of extragalactic 
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radiation are sufficiently well known to justify the assumption of gas neutrality for 
N(R I) > 2x 10 20 cm -2 . The close agreement between the mass per unit comoving 
volume of neutral gas in damped Lya systems and visible matter in current 
galaxies indicates that damped Lya systems comprise a significant neutral-gas 
reservoir for star formation at high redshift. 

(2) The comoving density of the neutral gas, £l g (z), declines by a factor of 
two between z ~ 3.5 and z ~ 2.3. While the evolution at < z < 2.3 is more 
uncertain, Q g (z) at z ~ 3.5 is a factor of three higher than at z = 0. 

(3) Damped Lya systems are metal-poor at all redshifts (see Table 2), but 
exhibit a metallicity "floor", [M/H] > —2.6 (Table 3), indicating a different en- 
richment history than that of the Lya forest. 

(4) The cosmic metallicity doubles every Gyr at z > 2, but the median [M/H] 
is sub-solar at z < 1.6. 

(5) From the large [Zn/Cr] ratios and the increase of the [Zn/Fe] and [Si/Fe] 
ratios with increasing metallicity we know that damped Lya systems exhibit 
evidence for depletion by dust and that the dust content is far lower than in the 
Galaxy. 

(6) The presence of a plateau in the [N/a] versus [a/H] plane near [N/a] ~ —0.7 
indicates a minimum age of 0.25 Gyr for damped Lya systems, which suggests 
they are not transient objects but instead probably have ages comparable to the 
Hubble time at the absorption epoch. 

(7) Ionized gas in damped Lya systems exhibits a different velocity structure 
than the neutral gas, unlike the agreement between the velocity structures of 
these two phases in the Galaxy. 

(8) H2 and other molecules are rarely present in damped Lya systems. Stud- 
ies of those systems exhibiting H2 absorption indicate the presence of an FUV 
radiation field with J v « 10~ 19 ergs cm~ 2 s _1 Hz -1 sr _1 , which resembles (a) 
the interstellar radiation field in the Galaxy and (b) J u predicted by the C II* 
technique. 

(9) The frequency distribution of the absorption velocity intervals, Av, has 
a median of 90 km s _1 . This property cannot be reproduced by single-disk 
CDM scenarios proposed so far, and is difficult to reproduce for sightlines passing 
through dwarf galaxies. Damped Lya systems with large values of Av exhibit a 
systematic absence of low values of [M/H] and high values of N(H I). 

(10) We cannot rule out the hypothesis that galaxies identified with damped 
Lya systems at z < 1.6 are drawn from a cross-section weighted sample of normal 
galaxies; i.e., an inflated populations of dwarfs is not required. 

(11) C II* A 1335.7 absorption is detected in about half of randomly-selected 
samples of damped Lya systems. The inferred [C II] 158 /im cooling rates indi- 
cate heating rates far in excess of those supplied by FUV background radiation, 
requiring a local heat source. The evidence accumulated so far suggests that the 
likely site of C II* absorption is CNM gas. 

Next, we describe critical unsolved problems in damped Lya research. 

(I) What is the median mass, M me ^, of the dark-matter halos containing damped 
Lya systems ? This is the critical diagnostic for discriminating among most hierar- 
chical models, in which M me( i < 10 9 M@, from hierarchical models with feedback 
or passive evolution models, in which M mec i > 10 11 M & . 

(II) Does the damped Lya luminosity function overlap that of Lyman Break 
galaxies? Partial Overlap is suggested by the luminosities of the few objects 
detected in emission. 
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(III) What are the properties of the interstellar gas in damped Lya systems? 
These are crucial for understanding whether or not the gas in which C II* ab- 
sorption is detected can support a CNM. 

(IV) Are stars forming in damped Lya systems when they are detected? 

(V) How are the star formation and accretion histories of damped Lya systems 
related? The C II* technique indicates that star formation depletes the neutral 
gas reservoir of damped Lya systems more rapidly than indicated by the decrease 
of £l g (z) with time at z > 2.3. Does this require accretion of neutral gas onto 
damped Lya systems at rates comparable to the star formation rates? 

(VI) What is the solution to the "missing metals" problem? Evidence for metal- 
enriched gas ejected from damped Lya systems, or for light emitted from compact 
"bulge" regions would help in deciding between these hypotheses. 

(VII) What is the intrinsic nucleo synthetic [Si/Fe] ratio in damped Lya sys- 
tems? Are the intrinsic abundances of damped Lya systems a-enhanced? 

(VIII) What is the cosmic metallicity of low-z damped Lya systems? Is the 
column-density weighted mean metallicity of low-z damped Lya systems biased 
by undersampling and by obscuration? 

(IX) How can we improve numerical simulations of damped Lya system evo- 
lution? The next steps involve more accurate modeling of star formation and 
mechanical feedback. 

A major goal of damped Lya system research is to give a clear and decisive 
answer to the question, "What is a damped Lya system?" Obviously this has not 
yet been accomplished. Rather, what we have found is that a significant fraction 
of damped Lya systems are a population of H I layers exhibiting many of the 
complexities of the ISM of the Galaxy. They clearly play an important role in 
the formation of galaxies and undoubtedly interact with other structures in the 
high redshift Universe through a variety of feedback mechanisms. Observations 
of damped Lya systems provide an amazingly rich data set that gives information 
about galaxy formation unavailable by other means. Specifically, observations of 
damped Lya systems are the only way to study in detail the neutral gas that 
gave rise to galaxies at high redshifts. We hope that the interplay between new 
observations and improved theoretical modeling will lead to significant insights 
into the process of galaxy formation. 
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DLA 




z DLA a 






kpc 


10- 17 (cg B ) 


SFR 
Diagnostic^ 


SFR Rcf 9 
M yr _1 


0458- 
0953+ 
2206- 


02 
-47A 
19A 


2.0395 
3.407 
1.9205 


2.0396 
3.415 
1.9229 


0.3±0.3 

<0.5 
1 


2.5±2.5 
<3.7 
S.4 


5 4 U ' 2 

0.8 
0.7±0.2 
26±3.0 


Lya 
Ly«,C 
C 


>1.5 1 
0.8^7.0 2 
26-*50 3 



1210+17 1.8918 0.25 2.1 <2.5 Ha <5.0 4 

1244-34 1.8590 0.16-0.24 1.4-2.0 <0.8 Ha <1.6 5 

8 DLAs 2.095-2.615 ssl.5 SBlO.O <9.0 Ha <30 6 



Table 1: Measurements of Emission from Damped Lya Systems 

a Absorption rcdshift of DLA 
fa Lya emission rcdshift 

c Displaccmcnt angle of emitter (or candidate) from QSO 
^Displacement distance of emitter (or candidate) from QSO 

e Line flux in units of ergs cm~^ s~ ^. First 3 entries are Lya and second 3 arc Ha 

^Diagnostic for obtaining SFR:Lya is Lya emission, C is FUV continuum, and Ha is from Ha emission. 

9 Reference:(l) Moller et al. (2004);(2)A. Bunker (2004), priv. comm; (3) Mollcr et al. (2002): (4) Kulkarni ct al. (2001);(5) Kulka- 
rni ct al. (2000);(6) Bunker ct al. (1999) 



Redshift Interval 
Property 0.0-1.6 1.6-4.5 
dM/dX .... 0.077±0.0lF 

W 3 n g {z) 0.96±0.23 a 0.92±0.21 b 
logi p'*(^) < -1-40 -0.70±0.28 c 
< Z > -0.81±0.032 -1.33±0.09 



Table 2: Global Properties 
a Due to large uncertainties in individual measurements, error in mean given by 
propagation of experimental errors. 

b Because of systematic decrease in dAf/dX and Q g (z) with time, error in mean 
determined by standard deviation. 

c Computed by averaging over the "WD low" and other dust models discussed in 
Wolfe et al. (2003a) (see their Table 1), and by assuming that both DLAs with 
detected and undetected C IF absorption have same mean SFR per unit area. 
If the SFR per unit area of the non-detections were significantly lower, p* would 
decrease by about 0.3 dex. 
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Property 


z 


N a 


x b 


%"med 


o A 


Min e 


Max-f 


Iog 10 JV(HI) 


z > 1.6 


199 


20.83 


20.60 


0.33 


20.30 


21.70 


[M/H] 


0.3 < z < 4.9 


130 


-1.11 


-1.48 


0.55 


-2.65 


0.04 


[Zn/Fe] 


0.7 < z < 3.3 


38 


0.54 


0.42 


0.25 


-0.01 


1.05 


[a/Fe] 


0.8 < z < 4.7 


70 


0.42 


0.38 


0.18 


0.03 


1.00 


Av,„ s 


1.7 < z 


95 


114. 


90. 


83.7 


16. 


430. 




1.7 < z 


75 


209. 


190. 


113.4 


20. 


528. 


logio/(H 2 ) 


2.0 < z < 3.4 


33 


-2.22 


< -5.93 


0.82 


< -6.98 


-0.64 


logio4 l 


1.7 < z < 4.2 


57 


-26.57 


< -26.93 


0.49 


< -27.69 


-25.35 




1.7 < z < 4.5 


39 


9.6 


5.4 


6.5 


<0.24 


23 


logi O '0* fe 


1.7 < z < 4.2 


40 


-1.95 


-2.20 


0.30 


< - 3.55 


-1.55 



Table 3: DLA Individual Summary 
"Number of DLAs in sample 

h Mean value. Upper limits excluded in computation. 
c Median value. Upper limits included in computation. 
^Sample dispersion 
e Minimum value. 

aximum value. 
9 Low-ion absorption velocity interval (in km s _1 ). 
''High- ion (C IV) absorption velocity interval (in km s _1 ). 
1 l c is in units of ergs s _1 H _1 

3 J u in units of 10 _19 ergs cm~ 2 s _1 Hz -1 sr _1 . Computed for WD low model in Wolfe 
et al. (2004) 

fe SFR per unit area for uniform disk model (in M Q yr _1 kpc~ 2 ). Computed for WD low 
model in Wolfe et al. (2004) 
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Figure 1: Keck/ESI spectrum of QSO PSS0209+0517 showing the Lya forest, a 
pair of damped Lya systems, and a series of metal-lines. The schematic labeling in 
the figure identifies several key features for the damped Lya system at z = 3.864. 
The absorption trough at A =5674 A corresponds to the damped Lya line at 
z = 3.667. 
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Figure 2: Example Voigt profile fits to two damped Lya systems of the sample 
from Prochaska et al. (2003b). The vertical dashed line indicates the line centroid 
determined from metal- line transitions identified outside the Lya forest. The 
dotted line traces the continuum of the QSO and the green and red lines trace 
the Voigt profile solution and the fits corresponding to lcr changes to iV(HI). The 
fluctuations at the bottom of the damped Lya absorption troughs indicate the 
level of sky noise. 
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Figure 3: The iV(H I) frequency distribution f(N,X) determined by Prochaska 
et al. (2005) for all damped Lya systems in the SDSS DR3-DR4 sample. Over- 
plotted on the data poinst are a single power-law, T function, and a double 
power-law. On the latter two are acceptable fits to the data. Plot taken from 
Prochaska et al. (2005). 
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Figure 4: Incidence of damped Lya systems per unit cosmological distance 
dM jdX (denote by £dla in the figure) as a function of redshift. The three 
data points at z = are all local measurements from 21 cm observations (Ryan- 
Weber, Webster, & Staveley-Smith 2003; Zwaan et al. 2001; Rosenberg & Schnei- 
der 2002). The curve overplotted on the data traces the evaluation of dN jdX in 
a series of 0.5 Gyr intervals. Plot taken from Prochaska et al. (2005). 
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Figure 5: Neutral gas mass density versus z from Prochaska et al. (2005). H 

I data at (a) z > 2.2 from SDSS DR3-4 survey, (b) < z < 1.6 from the Mg 

II survey of Rao, Turnshek, & Nestor (2004, priv. comm.), and (c) at z = 
(diamond) from Fukugita et al. (1998). Stellar mass density at z = (star) from 
Cole et al. (2001) and stellar mass density of Irr galaxies (+ sign) from Fukugita 
et al. (1998). Theoretical curves from Cen et al. (2003) (green), Somerville et al. 
(2001) (red), and Nagamine et al. (2004a) (blue: dotted is D5 model and solid is 
Q5 model) 
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Figure 6: Velocity profiles of metal-line transitions for (a) the metal-poor damped 
Lya system at z = 3.608 toward Q1108-07 and (b) the metal-strong damped Lya 
system at z = 2.626 toward Q0812+32. Grey indicates the range of flux between 
and 1. Red lines are blends due to other transitions. 
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Figure 7: Current summary of the metallicity measurements of the damped Lya 
systems as summarized in Prochaska et al. (2003a), Kulkarni et al. (2005), and 
Rao et al. (2005). The upper panel plots metallicities against redshift and the 
binned points indicate the cosmological mean metallicity with 95% c.l. uncer- 
tainty. The lower panel plots the metallicity versus look-back time. The over- 
whelming majority of observations are from t > 10 Gyr. 
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Figure 8: Plot of (a) gas-phase [Zn/Fe] observations against Zn metallicity for 
all damped Lya systems with high signal-to-noise echelle observations and (b) 
gas-phase [Si/Fe] observations against Si metallicity for all damped Lya systems 
with high signal-to-noise echelle observations. Damped Lya systems with upper 
limits to Zn were suppressed from panel (a). 
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Figure 9: Abundance pattern for the metal-strong damped Lya system at 
z = 2.626 toward Q0812+32. Because of the high metal abundance, [O/H] = 
—0.44, a dust correction is necessary, and in this case a conservative 'warm halo' 
correction (Savage & Sembach 1996) was applied. The dotted line traces the 
Solar abundance pattern scaled to match the oxygen abundance of the damped 
Lya system. 
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Figure 10: Histogram of the low-ion velocity widths for the current sample of 
damped Lya systems with HIRES, UVES, or ESI observations. The median 
Av value is 90km s" 1 and the distribution shows a significant tail to beyond 
200km s' 1 . 
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Figure 11: Comparison of the C IV velocity width with the low-ion velocity width 
for the current sample of damped Lya systems with HIRES or ESI observations. 
With only one or two exceptions, Aw civ >~ ^Low 
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Figure 12: £ c versus iV(H I) for sample of 52 damped Lya systems. Red data 
points are positive detections, blue are 2a upper limits, and green are 2a lower 
limits. Small stars depict positive detections from sightlines in the Galaxy ISM. 
Large star depicts [C II] 158 fim emission rate per H atom averaged over the disk 
of the Galaxy. The latter is about 30 times higher than the average of the DLA 
detections. 



